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ABSTRACT 



We explore the nature of variations in dust emission within an individual 
^ ■ galaxy using 3.6 - 160 /im Spitzer Space Telescope observations and 450 and 

850 ixm James Clerk Maxwell Telescope observations of the edge-on Sd spiral 
• galaxy NGC 4631 with the goals of understanding the relation between poly- 

^ ■ cyclic aromatic hydrocarbons (PAH) and dust emission, studying the variations 

■ in the colors of the dust emission, and searching for possible excess submillimeter 

^ ■ emission compared to what is anticipated based on dust models applied to the 



O ■ mid- and far-infrared data. PAH emission at 8 /im is found to correlate best with 

• hot dust emission at 24 fim on kiloparsec scales, although the relation breaks 

down on scales equal to hundreds of parsecs, possibly because of differences in 
2 ■ the mean free paths between the photons that excite the PAHs and heat the 
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dust and possibly because the PAHs are destroyed by the hard radiation fields in 
the centers of some star formation regions. The ratio of 8 fim PAH emission to 
160 fim cool dust emission appears to be a function of radius. The 70/160 and 
160/450 fim flux density ratios are remarkably constant even though the surface 
brightness varies by a factor of 25 in each wave band, which suggests that the 
emission is from dust heated by a nearly-uniform radiation field. Globally, we 
find an excess of 850 /im emission relative to what would be predicted by dust 
models. The 850 /im excess is highest in regions with low 160 /im surface bright- 
ness, although the strength and statistical significance of this result depends on 
the model fit to the data. We rule out variable emissivity functions or ~ 4 K 
dust as the possible origins of this 850 fim emission, but we do discuss the other 
possible mechanisms that could produce the emission. 

Subject headings: galaxies: ISM, infrared: galaxies, galaxies: individual (NGC 
4631) 

1. Introduction 

Before the launch of the Spitzer Space Telescope (Werner et al. 2004), research into the 
1 - 1000 yUm spectral energy distributions (SEDs) of nearby galaxies was greatly restricted 
by the angular resolution and signal-to-noise levels of far-infrared data. Most of the Infrared 
Astronomical Satellite (IRAS) and Infrared Space Observatory (ISO) data were only usable 
as global fiux density measurements. In the few cases where galaxies were resolved at far- 
infrared wavelengths, the spatial information was limited. With Spitzer data, it is now 
possible to examine 1 - 1000 fim SEDs for kiloparsec-sized regions in galaxies at distances 
up to 10 Mpc. These higher-resolution, higher-sensitivity data allow for new observational 
investigations into the properties of the polycyclic aromatic hydrocarbons (PAHs) and other 
dust across this large wavelength regime. In this paper, we will use Spitzer data to address 
two major issues related to dust emission that were not completely resolved using IRAS and 
ISO data. 

The first major question that can be asked is how best to describe the SEDs of discrete 
regions within individual nearby spiral galaxies, particularly in the context of describing 
the dust emission between 60 and 2000 /xm. Studies combining IRAS and ISO results with 
submillimeter or millimeter data suggested that a range of simple one or two-component 
blackbodies modified with different emissivity laws could describe the data. The exact 
descriptions, however, were confusing and contradictory. The major point of debate was 
how to describe the emission at submillimeter or millimeter wavelengths that exceeded what 
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is expected when extrapolating from the 60 - 200//m regime using one or a series of ~ 20 - 
30 K blackbodies modified by emissivity functions. Some studies have suggested that 
the power law describing the emissivity function varies among nearby galaxies (e.g Dunne 
et al. 2000; Bendo ct al. 2003). Other studies found that the submillimeter or millimeter 
excess could be described by a cold blackbody modified by a A^^ emissivity law. This cold 
dust temperature component could have a temperature as low as < 10 K (e.g Kriigel et al. 
1998; Siebenmorgen et al. 1999; GaUiano et al. 2003, 2005). 

Regan et al. (2004) were among the first to examine this issue using Spitzer data. They 
examined the SEDs of the central starburst ring and outer disk of NGC 7331 and found 
the 70 - 850 /xm dust emission from the ring was consistent with dust in the ~ 20 - 25 K 
range with an emissivity law of A"^. However, a limitation of the study by Regan et al. was 
the lack of significant submillimeter emission at a high signal-to-noise level outside of the 
ring. They therefore could not determine whether the description of the dust emission in 
the center of NGC 7331 apphes throughout the galaxy's disk. 

This issue is particularly important in terms of modeling the dust emission and esti- 
mating the dust mass within nearby galaxies. If dust as cold as ~ 6 K is present, it will 
only be a significant source of emission at submillimeter wavelengths, but it may contribute 

substantially to the mass of the interstellar medium. If the dust emissivity is variable then 
dust models need to be adjusted to account for this variability. Without these results, the 
full picture of dust emission in nearby galaxies is incomplete. 

The other major question is to determine how PAH emission, particularly the 7.7 fj,m 
feature, is correlated with dust emission on scales of hundreds of parsecs. Results from 
ISO gave differing results. Forstcr Schreibcr ct al. (2004), among others, found a strong 
correlation between 7 and 15 fira emission in galaxy disks, thus demonstrating a correla- 
tion between PAH and hot (~ 100 K) dust emission (which is commonly associated with 
transiently-heated very small grain emission but which may also represent emission from 
grains that are in thermal equilibrium at ~ 100 K). Haas et al. (2002), however, found that 
8 iim PAH emission was correlated more closely with 850 iim emission from cool (~ 25 K) 
grains than with the 15 /xm emission from hot dust. Using 5pzteer observations of NGG 300, 
HcIoTi ct al. (2004) found that 8 /im emission originated from the rims of star formation re- 
gions within the galaxy whereas the 24 /xm emission was more strongly peaked in the centers 
of the star formation region, which suggested a correlation on kiloparsec scales but not on 
smaller scales. Dale et al. (2005), however, found that 8 fim emission from kiloparsec sized 
regions within M 81, M 51, and NGC 7331 appeared to be more closely correlated with the 
total infrared flux than with the 24 /im flux density. In a separate analysis of the Spitzer 
data for M 81, Perez- Gonzalez et al. (2006, in press) also found that the 8 /xm luminosity 
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was closely correlated with the total infrared luminosity on kiloparsec scales, although they 
did not examine the relation between 8 and 24 /im emission. 

Identifying how PAH emission is correlated with dust emission will lead to a better 
understanding of how it can be used as a star formation indicator, as has been suggested 
by the results of Roussel et al. (2001) and Forster Schreiber et al. (2004) (although the 
contradictory results of Calzctti et al. (2005) should also be noted). If PAH emission is 
correlated with 24 fim hot dust emission on scales of hundreds of parsecs or kiloparsecs, then 
it should be as effective a star formation indicator as 24 fim emission on those spatial scales. 
However, if PAH emission is more closely correlated with 160 fim cool dust emission, which 
may trace extended cirrus emission unrelated to star formation (e.g. Hinz et al. 2004), then 
the PAH emission should be used with greater caution as a star formation tracer. 

Among the data available in the Spitzer Infrared Nearby Galaxy Survey (SINGS; Ken- 
nicutt et al. 2003), the data for NGC 4631 are optimal for studying the complete SED of 
the dust emission from mid-infrared to submiUimeter wavelengths. NGC 4631 is a nearly 
edge-on (inclination 85 deg; TuUy 1988) Sd galaxy at a distance of 9.0 Mpc with an optical 
disk of 15'. 5 X 2'. 7 (de Vaucouleurs et al. 1991). The edge-on nature of the galaxy makes 
it a high surface brightness source that is easier to detect at all wavelengths, particularly 
submillimctcr wavelengths. However, the edge-on orientation adds a geometry effect to the 
observed surface brightness variations; regions may appear bright because they are intrin- 
sicly luminous or because they represent regions with high column densities. Nonetheless, in 
terms of this study, the advantages of the edge-on orientation outweigh the disadvantages. 

In addition to the Infrared Array Camera (IRAC; Fazio ct al. 2004) and Multiband 
Imaging Photometer for Spitzer (MIPS; Riekc et al. 2004) data taken as part of the SINGS 
legacy program, deep archival SCUBA (Holland et al. 1999) 450 and 850 fim data covering 
most of the optical disk (previously published in Alton et al. 1999; Stevens ct al. 2005) are 
available for this galaxy. These data almost completely cover the extended dust emission in 
this galaxy and allow for extraction of 3.6-850 /im SEDs from regions throughout the disk. 
Mid-infrared-to-submillimeter SEDs of this galaxy have been studied before (Bendo et al. 
2003; Dumke et al. 2004; Stevens et al. 2005). However, these previous studies have rehed on 
IRAS and ISO data that have lacked the resolution, the sensitivity, or the spatial coverage to 
study anything more than the global SED or the SED of the center of the object. With these 
new data, we can study infrared color variations across the entire optical disk of the galaxy 
and examine the SEDs within discrete regions throughout the disk. These results will lead 
to answers as to how to characterize far-infrared to submiUimeter dust emission throughout 
a spiral galaxy and how PAH and dust emission at different wavelengths are interrelated. 

In Section 2, we present information on the observations and data reduction for this 
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galaxy. We briefly discuss the images in Section 3. In Section 4, we present simple color 
information from flux densities integrated over set apertures distributed throughout the plane 
of the galaxy. Then, in Section 5, we present the 3.6 - 850 /xm spectral energy distributions 
of the galaxy as a whole and of discrete regions within the galaxy. We then summarize the 
results in Section 6. 



2. Observations and Data Reduction 

2.1. IRAC (3.6, 4.5, 5.7, and 8.0 //m) Data 

The 3.6 - 8.0 //.m data were taken with IRAC on the Spitzer Space Telescope on 2004 
May 23 and 26 in IRAC campaign 5. The observations consist of a series of dithered 5' x 5' 
individual frames that are offset 2. '5 from each other. The two separate sets of observations 
allow asteroids to be recognized and provide observations at orientations up to a few degrees 
apart to ease removal of detector artifacts. Points in the center are imaged eight times in 30 
s exposures. The full-width half-maxima (FWHM) of the point spread functions (PSFs) are 
l."6 — 2."0 for the four wavebands. 

The data are processed using a special SINGS IRAC pipeline. First, a geometric distor- 
tion correction is applied to the individual frames. Data from the second set of observations 
arc rotated to the same orientation as the first set of observations. Bias structure is sub- 
tracted from the data by subtracting a bias frame from each frame. This bias frame is made 
by median combining all data frames for the observations; the median filtering effectively 
removes any target signal from the data, leaving only the bias signal. Next, the image off- 
sets are determined through image cross-correlation. Following this, bias drift is removed 
by determining the relative offset between overlapping regions of individual frames and sub- 
tracting the offset. Finally, cosmic ray masks are created using standard drizzle methods, 
and final image mosaics are created using a drizzle technique. The final pixel scales are set 
at ~ 0".75. A final residual background is measured in small regions outside the target that 
are free of bright foreground/background sources, and then this residual is subtracted from 
the data. The uncertainty in the calibration factor applied to the final mosaics is 10%. 

2.2. MIPS (24, 70, and 160 /xm) Data 

The 24, 70, and 160 /xm data were taken with MIPS on the Spitzer Space Telescope on 
2005 January 24 and 26 in MIPS campaign 18. The observations were obtained using the 
scan-mapping mode in two separate visits to the galaxy, which has the same benefits as for 
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the IRAC data. As a result of redundancy inherent in the scan-mapping mode, each pixel 
in the core map area was effectively observed 40, 20, and 4 times at 24, 70, and 160 /xm, 
respectively, resulting in integration times per pixel of 160 s, 80 s, and 16 s, respectively. 
The FWHM of the PSFs are 6// at 24 /xm, 18// at 70 /xm, and 40// at 160 /xm. 

The MIPS data were processed using the MIPS Data Analysis Tools (Gordon et al. 2005) 
version 2.96. Since the processing for the 24 ^m data differed from the 70 and 160 /xm data, 
the 24 /im data processing will be discussed separately from the 70 and 160 /im wavebands. 

First, the 24 /im images were processed through a droop correction (that removes an 
excess signal in each pixel that is proportional to the signal in the entire array) and were 
corrected for non-linearity in the ramps. The dark current was then subtracted. Next, scan- 
mirror-position dependent fiats (created from all SINGS MIPS off-target data from MIPS 
campaign 18) were applied to the data. Following this, scan- mirror-position independent flats 
(made from off-target data frames within each visit's scan map) were applied. A correction 
to the offsets between the four channels reading out the 24 //m detector data was applied to 
each data frame. Latent images from bright sources, erroneously high or low pixel values, 
and unusually noisy frames were also masked out. For background subtraction, the median 
background value as a function of time was determined for each scan leg and then subtracted. 
The data from both visits were then mosaicked together, and a final residual background 
offset was measured outside the optical disk of the galaxy and subtracted. We applied the 
SIO calibration factors in the MIPS Data Handbook^^ {Spitzer Science Center 2006) to the 
final mosaics. The calibration factor has an uncertainty of ~10%. Color corrections, which 
are smaller than the calibration uncertainties, are not applied to these data except for making 
comparisons to IRAS 25 /xm data. 

For the 70 and 160 fim data, ramps were fit to the 70 and 160 fim reads to derive 
slopes. In this step, readout jumps and cosmic ray hits were also removed, and an electronic 
nonlinearity correction was applied. Next, the stim fiash frames taken by the instrument 
were used as responsivity corrections. The dark current was subtracted from the data, an 
illumination correction was applied, and then short term variations in drift were subtracted 
from the data (which also subtracted the background). Following this, erroneously high or 
low pixel values were removed from the data. The data from both visits were then mosaicked 
together, and a residual offset measured in regions staggering the target was subtracted from 
the final maps. We applied the SIO calibration factors in the MIPS Data Handbook, which 
have uncertainties of ~20%, to the final mosaics. These uncertainties can also be treated 
as encompassing the uncertainties from nonlinearity effects in the 70 /^.m wave band. The 
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color corrections are smaller than the calibration uncertainties and are not applied to these 
data. Note that the 70 /im data were affected by a latent image effect that created negative 
streaks to the south of the central ~ 2' of the data. However, we will not be using this region 
in the analysis. 

To check the accuracy of the flux density measurement at 24 /im, we can compare the 

integrated flux density with the 25 //m integrated flux density measured from IRAS data. The 
24 fim global flux density measured here is 8.0 ±0.8 Jy. The MIPS Data Handbook suggests 
a color correction term of ~ 0.95 for blackbody emission from a 100-200 K somcc{Spitzer 
Science Center 2006), which is what is anticipated for the dust emission that dominates the 
24 yum wave band. Dividing the observed flux density by this factor gives a color-corrected 
24 iim flux density of 8.4 Jy. Rice et al. (1988) published a measured 25 iim flux density 
of 9.65 ± 1.4 Jy and a color-corrected flux density of 8.35 Jy. Given the uncertainties in 
the MIPS and IRAS measurements, the color-corrected flux densities are in good agreement 
with each other. 



2.3. SCUBA (450 and 850 /im) Data 

The 450 and 850 iim data were taken from the JCMT SCUBA archives maintained by 
the Canadian Astronomy Data Centre. The data used were taken in jiggle map mode. Each 
jiggle map consists of data taken in a ~ 2.3' region. During the observations, the secondary 
shifts through a 64-point pattern so that the bolometers may fully sample the image plane. 
Additionally, the secondary chops off-source so as to measure the background during the 
observations. In the data used here, the chop was 120'-180' in the north-south direction. 
The FWHM of the PSFs are 8// and 14// at 450 and 850 /^m, respectively. 

Data from projects with identification numbers U43 (5-6 Oct 1997), M98BU19 (8-9 
Mar 1998), M99AI16 (14 May 1999), and M00BU51 (14 Jan 2001) were used to build the 
submillimcter images. These data consist of simultaneous observations at 450 and 850 /xm, 
although the data from 14 May 1999 and one frame of data from 14 Jan 2001 were not 
usable for constructing the 450 /xm map. Note that the SCUBA narrow-band filters were 
used for the 1997-1999 observations and the wide-band filters were used for 2001 data; except 
for some improvements in the signal-to-noise in the 450 /im wide-band filter data, the two 
filter sets are otherwise functionally similar and treated that way in the data processing. All 
together, the data consist of 182 integrations (194 min) of on-target observations (of which, 
149 integrations (159 min) were usable 450 /xm data). 

The data were processed with the SCUBA User Reduction Facility (Jenness & Light- 



-8- 



foot 1998). The data were first flatfielded and corrected for atmospheric extinction. Noisy 
bolometers were removed next, followed by spike removal. Then the background signal was 
subtracted using the signal from several bolometers at the north and south edges of each 
frame that were selected so as to avoid noisy pixels and any possible negative images (which 
are an artifact related to the chopping performed during the observations). The data were 
calibrated using observations of the planets Mars or Uranus or the submillimeter standard 
CRL 618 and then regridded onto the sky plane. An additional background subtraction was 
then performed to remove any residual offset. A 6" border at the edge of the observed field 
was masked out so as to remove anomalously "hot" or "cold" pixels that are artifacts of 
the sky subtraction. These anomalous pixels are identifiable by the amplitude of the signal, 
which is greater than the signal from the target, and by the shapes of their radial profiles, 
which are narrower than the PSFs of true sources. Additional anomalous pixels were iden- 
tified by eye and masked out. Negative pixels near the south edge of the map just cast of 
the nucleus were identified as a negative image of the plane of the galaxy and masked out in 
the final map, and the region around a noisy bolometer at the west end of the 850 fim map 
was masked out as well. Typical uncertainties in the calibration are ~ 25% at 450 /im and 
~ 10% at 850 //m. 

To test the data reduction, the flux densities can be compared to the 450 and 850 /im 
flux densities measured in the same data by Stevens et al. (2005) and the 870 /xm flux density 
measured by Dumke et al. (2004) with the 19-channel bolometer array at the Heinrich-Hertz 
Telescope (HHT). The flux densities measured with the archival SCUBA data in this paper 
(without adjustments for CO, free- free, or synchrotron emission) are 24±6 Jy at 450 /xm and 
4.9 ±0.5 Jy at 850 fim. Stevens et al. (2005) measured statistically identical 450 and 850 
global flux densities (without any adjustments) from the same data, thus demonstrating that 
the data reduction techniques used here can accurately reproduce the data. Dumke et al. 
(2004) measured the 870 fim fiux density (without any adjustments) to be 3.8 ±0.6 Jy. This 
is only a ~ 2a discrepancy between the SCUBA and HHT measurements. The SCUBA data, 
however, have a higher signal-to-noise ratio throughout the disk than the HHT data. The 
plane of the disk from ~ 200" east of the nucleus to ~ 300" west of the nucleus is detected at 
the 3(7 level in the unsmoothed 850 /im map presented here. In the smoothed 870 /im data of 
Dumke et al., however, some regions inside the disk of the galaxy with signal-to-noise levels 
of less than 2a are visible, particularly to the east of the center. In the central 135", where 
the signal-to-noise levels of both the SCUBA data and the Dumke et al. data are greater 
than 5o", the two maps give results that agree with each other. We measure the fiux density 
to be 1.8 ± 0.2 Jy within this region, whereas Dumke et al. obtained 1.8 Jy. Therefore, 
we conclude that the disagreement is for the flux density measurements in the outer disk. 
Because the signal-to-noise levels of the SCUBA data are superior to the Dumke et al. data. 
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we will rely only on the SCUBA data for 850 //m flux densities. (The signal-to-noise ratio 
of the 1230 //m data in Dumke et al. are higher than 3ct throughout the disk, so these data 
should be more reliable. We will use these 1230 /im data in constructing global SEDs.) 

2.4. Convolution of the Data 

To compare data between wavelengths of different resolutions on scales smaller than an 
arcminute, it is best to convolve the data with kernels that will match not only the FWHM of 

the PSFs but also the shapes of the PSFs. Wc arc using convolution kernels that will convert 
an input PSF into a lower resolution output PSF using the ratio of Fourier transforms of 
the output to input PSFs. High frequency noise in the input PSF is suppressed when these 
kernels are created. See K. D. Gordon et al. (2006, in preparation) for details. 

For comparing flux densities between two wave bands, we match the resolutions of the 
data to either the 18// resolution of the 70 fim data or the 40// resolution of the 160 /xm data. 
In a comparison of PAH 8 fim and 24 fim maps in Section 4.2, we degraded the IRAC data 
to the 6// resolution of the 24 /xm data. Note that the kernel for degrading the resolution of 
the 450 /xm maps to 18// was not available at the time of this writing. 

3. Images 

The global 3.6-850 images of these data are presented in Figure 1. The PAH and 
dust emission structures from 5.7 fim to 850 fim all appear very similar, and they do trace 
similar structures seen in the 3.6 and 4.5 /xm bands, which contains the Rayleigh- Jeans side 
of the blackbody starlight emission. Most of the emission lies on the major axis of the optical 
disk of the galaxy. No bulge is visible, even in the 3.6 and 4.5 /xm bands. 

The brightest 3.6-850 /im source is a central structure with a diameter of ~ 2' that looks 
like an edge-on ring. Magnified images of this structure are shown in Figure 2. Both the east 
and west lobes of this central structure as well as the nucleus are the three single brightest 
sources in the central region, although the relative brightness of the nucleus compared to 
the lobes varies. The nucleus is the brightest 3.6-24 //m source, but the east lobe in the 
central structure is brighter at 70-850 /im. Aside from this major variation in the colors, 
some other more subtle variations can be seen. Two knots of infrared emission to the west 
of the nucleus are stronger 24 and 70 fim sources than the filamentary structure just east of 
the nucleus, but this filament to the east of the nucleus is a brighter 450 and 850 fim. source 
than the knots to the west of the nucleus. Despite these minor color variations, however, all 
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wavebands effectively trace the same structure. 

Outside of the nucleus, all wavebands from 5.7 to 850 /xm still trace similar structures. A 
series of knots loosely connected by some filamentary dust structures can be seen throughout 
the disk of this galaxy. The appearance of the disk is notably asymmetric. The infrared 
emission clearly decreases more sharply to the east of the nucleus than it does to the west. 
The knots to the west of the central structure all lie in a single line. The dust structures 
to the east, however, are more jumbled in appearance. Some of the emission above and 
below the major axis of the optical disk of the galaxy may origniate from outside the plane 
of the galaxy or near the edge of the stellar disk, but it is difficult to distinguish which 
interpretation is correct from these data alone. 

The interacting elliptical galaxy NGC 4627 can be seen to the northwest of the center 
outside the plane of the galaxy. The infrared emission from this galaxy is mostly stellar; in 
these data, it is brightest at 3.6 /im, and it generally decreases in brightness towards longer 
wavelengths. Faint emission is still vaguely visible at 24 /xm, but none is visible at 160 /im. 
No significant dust structures connecting NGC 4627 and NGC 4631 are visible. 

4. Color Veiriations across the Dust Disk 

4.1. Measured Flux Densities within Discrete Regions 

To study the infrared color variations within NGC 4631, we have extracted flux densities 
from a series of circular regions within the optical disk of the galaxy. These regions were 
selected by eye to include all the bright knots in the 24-70 /im images as well as a few 
other prominent infrared structures. Two different diameters for the regions are used: 15" 
(to approximately match the FWHM of the 70 and 850 fim PSFs) and 40" (to match the 
FWHM of the 160 fim PSF). At a distance of 9.0 Mpc (Kennicutt et al. 2003), the 15" and 
40" angular sizes of these regions correspond to diameters of 650 and 1700 pc. In these 
comparisons, the data are flrst convolved with the kernels described in Section 2.4 so that 
all the image resolutions match. 

The regions in which the flux densities are measured are shown in Figure 3. Right 
ascensions, declinations, and flux densities are listed in Tables 1 and 2. Uncertainties in 
the calibrations are listed as percentages, and uncertainties from background noise (mea- 
sured in off-source regions of the sky) are fisted in Jy. Note that the uncertainties from 
the background noise is negligible compared to the calibration uncertainty except in the 
submillimeter bands. Also note that the uncertainties of the mean background levels (which 
are subtracted from each waveband) correspond to uncertainties in the flux densities that is 
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negligible compared to the uncertainties from background noise. For comparing the ratios 
of flux density measurements from two individual wavebands, only the background noise is 
relevant. In plots of flux density versus the ratio of flux densities, changes in the calibration 
factors will have systematic effects on individual wavebands that will only change the offsets 
in the plots. Since the goal is simply to identify whether the ratio varies as a function of 
flux density, any offset introduced by a calibration uncertainty does not change the scientific 
results. For composing SEDs in Section 5, however, the calibration uncertainties will also 
be used, since changes in calibration factors could change the shapes of the SEDs as well as 
the models fit to the data and the physical interpretation of the results. 

The IRAC measurements require the application of aperture corrections to account 
for a the scattering of hght through the detector material (Reach et al. 2005), a problem 
which, among the data in this paper, is unique to the IRAC data. Note that this IRAC 
aperture correction is not the same as an aperture correction that accounts for the fraction 
of the PSF created by the optics of the telescope that falls outside a measurement region 
(which docs not need to be applied if the data being compared are convolved to matching 
resolutions). Aperture corrections are not applied to the flux densities measured in the 15" 
regions because the flux densities measured in these regions are used for only studying color 
variations (so aperture corrections would only offset the data points in the plots), because no 
definitive aperture corrections have been published at the time of this writing, and because 
preliminary work on IRAC aperture corrections implies that the aperture corrections will 
be approximately equal to or less than the calibration uncertainties (T. Jarrett, private 
communication). However, since the flux densities measured in the 40" regions are used for 
modeling the SEDs (where precise calibration is more important), the "infinite" aperture 
corrections of Reach et al. (2005) have been applied to the flux densities measured in the 
40" regions. The flux densities listed in Table 2 include these corrections. 

To use the 8 //m band as a tracer of PAH emission, we need to subtract the stellar 
continuum emission. To do this, we will apply the formula 

f^,{PAH8fim) = f^{8fj.m{total)) - 0.232/^(3.6/xm) (1) 

from Helou et al. (2004) to the measurements in Tables 1 and 2. We use the term "PAH 
8 /um" to indicate where the continuum has been subtracted. Note that this continuum 
subtraction only removes the contribution of starlight. In regions with low PAH emission, 
a significant fraction of the non-stellar emission may originate from hot dust (see Cannon 
et al. 2006, , for example). Therefore, the actual variations PAH emission may be stronger 
than what is represented by the continuum-subtracted "PAH 8 //m" . 
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4.2. Compcirison between PAH and Hot Dust Emission 

First we will investigate how 8 /xm PAH emission and 24 /xm hot dust emission are 
interrelated. 

Figure 4 shows how PAH 8 fim and 24 /im flux densities measured within the 15" 
(650 pc) regions in Table 2 and the 40" (1700 pc) regions in Table 2 are related to each 
other. The best fitting functions are 

(-0.03 ± 0.04) + (0.032 ± 0M9)log{f24^cm/ Jy) (15"regions) (2) 

(-0.13 ± 0.03) + (0.08 ± 0.02)log{f2A;,m/Jy) (40"regions). (3) 

Note that the correspondence between the two wavebands is very good when the flux densities 
are measured within 40", but the relation begins to break down in 15". In both aperture 
sizes, the (PAH 8)/24 fim flux density ratio does not vary significantly with the 24 fim flux 
density. However, the standard deviation of the (PAH 8)/24 /xm flux density ratio is a factor 
of 2 higher in the 15" data compared to the 40" data, partly because more regions with higher 
(PAH 8)/24 iim flux density ratios are probed by the 15" regions. The map of the ratio of 
PAH 8 fjLm to 24 //m surface brightness in Figure 5 shows exactly why the data exhibit this 
scatter. The map reveals a significant amount of variation in the ratio within the disk of the 
galaxy on sub-kiloparsec scales. 24 ^m emission appears to peak more strongly than PAH 
emission in the centers of the dusty knots of star formation in this galaxy. The PAHs may 
be destroyed in the centers of these knots, whereas the hot dust emission rises because the 
dust is heated more strongly by the intense radiation fields in these knots. This phenomenon 
had also been reported for NGC 300 in Helou et al. (2004), and further Sj^zfeer observations 
may reveal that this phenomenon is common among nearby galaxies. 

To further examine the difference between the distribution of PAH and hot dust emis- 
sion, we extracted the PAH 8 yum and 24 yum radial profiles of multiple individual regions: the 
center of the galaxy, the lobes at either end of the central structure, and the seven brightest 
sources outside the center. These radial profiles are shown in Figure 6. As can be seen in 
these figures, the PAH 8 //m emission is clearly broader than the 24 //m emission. The PAH 
8 //m data used to create these radial profiles have been convolved to match the PSFs of 
the 24 /im data, but in several of the regions considered here, the PAH 8 //m radial profiles 
extracted from unconvolved data are still broader than the 24 /im radial profiles. 

In summary, these data demonstrate that, on kiloparsec scales, PAH and hot dust 
emission are correlated, but on smaller scales, the two can be spatially separated from each 
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other. The breakdown of the correlation on small scales can be physically interpreted in 
two ways. One possibility is that the primary source of heating for the dust that produces 
the 24 /xm emission may be ultraviolet light that is absorbed close to the centers of the star 
formation regions, but the PAHs can be excited by visible light (Li & Drainc 2002) that has 
a greater mean free path. The PAH emission would therefore appear to be more extended 
than the 24 /xm emission. The other possibility is that PAHs may be destroyed by the hard 
radiation in the centers of star formation regions (e.g. Madden et al. 2006), which is why 
the PAH 8 fjLm emission does not appear to be as peaked as the 24 //m emission. Either of 
these mechanisms would be an adequate explanation for the effects observed here, and both 
mechanisms may be involved in creating the observed color variations. 



4.3. Compcirison of PAH and Hot Dust Emission to Cool Dust Emission 

Next, we will investigate how the PAH 8 /im emission and 24 /im hot dust emission are 
related to cooler dust emission at 70 and 160 //m. 

Figure 7 shows comparisons of hot dust emission at 24 fim to the Wicn side of the cool 
dust emission at 70 /xm and the Rayleigh- Jeans side of the peak of the cool dust emission at 
160 iim. Both plots show a general trend where the 24/70 yum and 24/160 jim flux density 
ratios increase with surface brightness. The functions that best fit these data are 

log (Ip^] = (-1.007 ± 0.012) + (0.237 ± O.OU)log{f7o^m/ Jy) (4) 

V /70/im / 

and 

log (^^] = (-1.70 ± 0.02) + (0.24 ± 0.03)log{Ueo^.J Jy)- (5) 

Because of the edge-on nature of this galaxy, it is unclear as to whether this color variation is 
linked to variations in the dust heating with the intrinsic source luminosity (as anticipated, 
for example, by the models of (Dale et al. 2001) and (Li & Draine 2001)) or if the variation 
simply demonstrates that locations with significant hot dust emission are associated with 
large column densities of cool dust. However, if the color variation is finked to the cool dust 
column densities, then the amount of dust integrated over the line of sight is not physically 
linked to the hot dust, so some significant scatter should be visibile in the relations between 
surface brightness and color. The absence of such scatter in Figure 7 implies that the color 
variations are more closely related to dust heating, although column effects may still play a 
significant role. 

Note the two data points in the plot of the 24/70 iim fiux density ratio versus 70 iim 
fiux density that deviate from the general trend. The high surface brightness deviant data 
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point represents the nucleus, and the other deviant data point represents the point-hke 
source at right ascension 12:42:21.4 and dechnation +32:33:05 located above the major axis 
of the optical disk of the galaxy to the east of the center. These regions may deviate 
from the trend because the star formation in these regions is strong enough yet compact 
enough that the infrared colors deviate significantly from the rest of the disk, that something 
other than star formation (possibly an active galactic nucleus in the center, for example) 
is heating the dust, or that the surface brightness of stars is high enough that the stellar 
emission in the 24 /im band is non-negligible. Also note that the extraplanar region may be 
a foreground/background source not associated with NGC 4631. 

Figure 8 shows comparisons of PAH 8 /xm emission to the large grain emission at 160 /xm. 
For comparison with Haas et al. (2002), we also include a comparison of PAH 8 /im emission 
to 850 /im emission, although the 850 /im wave band may not trace the same dust as the 
160 fjim wave band (see Section 4.4). These data show that the (PAH 8)/160 //m and (PAH 
8)/850 fim flux density ratios increase with surface brightness, which means that the PAH 
emission is higher in systems with higher surface brightnesses. The functions that best fits 
these data are 



These trends clearly show that the PAH emission does not have a one-to-one correspondence 
with the 160 or 850 iim emission in this galaxy, which contradicts the results of Haas et al. 
(2002). Note that Haas et al. observed variations in the PAH 8/850 iim ratio of less than 
a factor of 2, whereas the ratios presented here vary by more than a factor of 10. At least 
on kiloparsec scales within this galaxy, PAH emission is more closely related to hot dust 
emission than to cool dust emission. 

The variations in the ratio of PAH emission to large grain emission could arise from either 
luminosity effects (where infrared-luminous regions produce relatively more PAH emission) 
or radial effects (where fewer PAHs are found at larger radii). To investigate this further, 
wc compared plots of the total infrared luminosity (calculated using the MIPS 24, 70 and 
160 /xm data and equation 4 in Dale & Helou 2002), the 24/160 /xm flux density ratio, and 
the (PAH 8)/160 jim flux density ratio as a function of radius for all regions listed in Table 2 
that lie along the major axis of optical disk of the galaxy. These radial plots are shown in 
Figure 9. The plots show that the 24/160 /xm and (PAH 8)/160 /xm flux density ratios both 
peak in the center, which is the most infrared-luminous region in the galaxy. Outside the 
center, however, a correspondence can be found between local peaks in the 24/160 /im flux 
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density ratio and local peaks in the total infrared emission, which indicates that variations 
in the 24/160 fim flux density ratio are linked to infrared surface brightness. In contrast, 
the (PAH 8)/160 /im flux density ratio monotonically decreases with radius. No local peaks 
in the (PAH 8)/160 fira flux density ratio corresponding to local peaks in the total infrared 
emission are found. The variation in the (PAH 8)/160 /xm flux density ratio is therefore a 
radial effect. 

At first, it appeared that the variations in the (PAH 8)/160 fim flux density ratio may 
be related to metallicity since metallicity also generally varies with radius. Variations in the 
strength of PAH emission with metallicity have been found previously in observations of the 
integrated SEDs of galaxies (Thuan et al. 1999; Dale et al. 2005; Engelbracht et al. 2005; 
Galliano et al. 2005). However, Otte et al. (2002) find variations in the nitrogen abundance 
in this galaxy that are smaller than the variations in the (PAH 8)/160 pLm flux density ratio 
found here, though the examined region does not span the entire optical disk of the galaxy. 
Moreover, based on the data of Engelbracht et al. (2005), the variation between PAH and 
dust emission is expected to look like a step function, not a smooth gradient, although the 
edge-on orientation could be responsible for smoothing out the emission so that this step 
function disappears. Therefore, we are hesitant to indicate that the radial variations in the 
(PAH 8)/160 iim flux density ratio are connected to metallicity. The cause of the radial 
variations needs to be examined further. 



4.4. Variations in Lcirge Grain Emission 

Color variations in the large grain emission from 70-850 /im show how the color tem- 
peratures of the large grains vary with surface brightness and may also reveal hints of where 
the coolest dust emission is located within the galaxy. 

Figure 10 shows how the 70 fira, 160 //m, and 450 //m emission are interrelated. The 
functions that best fit these data are 

log (p^) = (-0.444 ± 0.013) + (0.010 ± OM6)log{^eo^.m/Jy) (8) 

V JieOixm J 

and 

log (Im^ ^ (0.96 ± 0.03) + (-0.06 ± 0.05)Zo^(/450Mm/ Jy)- (9) 

V /450/im / 

Functionally, this is an indicator of the heating of the large dust grains. Virtually no variation 
in the 70/160 //m or 160/450 //m colors is apparent, even though the flux densities vary by 
more than an order of magnitude. At first, this appears to contradict the results of dust 
models such as the semi-empirical dust models of Dale et al. (2001) and the physical dust 
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models of Li & Draine (2001), which anticipate that these colors should vary dramatically 
with changes in the radiation fields or surface brightnesses of the regions. However, the 
absence of variations in the 70/160 fim. and 160/450 fim. flux density ratios could indicate 
that the interstellar radiation fleld that heats the large grains does not signiflcantly vary 
throughout the disk of this galaxy. The cause of the variations in surface brightness would 
therefore be related to variations in the mass of the dust producing the infrared radiation. 
However, since this galaxy is viewed edge-on, the color variations in the disk may be averaged 
out along the line of sight, although we may still expect some differences between infrared- 
bright and infrared-faint regions. 

For comparison, ISO observations of targets such as M 51, M 101 (Hippelein et al. 1996), 
NGC 6946 (Tuffs et al. 1996), and the disk of M 31 (Haas et al. 1998) show that far infrared 
colors do not vary with surface brightness. However, some Spitzer observations of nearby 
galaxies such as M 81 (Perez-Gonzalez et al. 2006, in preparation), M 51, NGC 7331 (Dale 
et al. 2005), and NGC 55 (Engelbracht et al. 2004) do show variations in the 70/160 fxm 
color with surface brightness. In light of these results, the invariance of the 70/160 //m colors 
found here might be unusual. More observations of far-infrared color variations within other 
galaxies are needed to interpret the galaxy-to-galaxy differences. 

Figure 11 shows how 850 fim emission varies with the large grain emission at 70 fim, 
160 /im, and 450 /im. The functions that best fit these data in logarithm space are given by 

(1.49 ± 0.11) + (0.21 ± Om)log{fs5o^.m/Jy): (10) 
(2.14 ± 0.14) + (0.47 ± O.U)logifs5o^m/Jy), (11) 

(1.06 ±0.08) + (0.37±0.07)/o5(/850Mm/^2/). (12) 

If the 160 /im, 450 /xm, and 850 fj,m emission sample the Rayleigh- Jeans side of ~ 25 K dust 
emission and if the emissivity law docs not vary, then the 160/850 fim and 450/850 fim flux 
density ratios should remain constant. Furthermore, since the 70 fim, 160 fira, and 450 //m 
surface brightnesses are well correlated, the dust heating should be uniform, and the ratio of 
850 /xm emission to emission at shorter wavelengths should not change signiflcantly. Despite 
these expectations, strong variations in 850 /xm relative to 70-450 /xm emission is visible in 
these data. As the surface brightness of the regions decreases, the 850 //m emission increases 
relative to the 70-450 iim emission. These trends suggest that the 850 iim emission cannot 
also be described as originating entirely from dust with a temperature of ~ 25 K and an 
emissivity of A~^. 
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4.5. Conclusions on Color Variations 

In summary, this comparison of flux densities measured in different wavebands has 
yielded interesting results. The PAH 8 /im and 24 /im hot dust emission are found to be 
correlated on kiloparsec scales, but the correlation breaks down on the sub-kiloparsec level. 
The 24 iim hot dust emission is found to increase relative to the 70 and 160 iim cool dust 
emission as the far-infrared surface brightness increases, as had been anticipated based on 
previous dust models. The PAH 8 //m emission is also found to increase relative to the 160 
and 850 /im wave bands as the far-infrared surface brightness increases, partly because of 
variations with radius. The 70/160 /xm and 160/450 /xm flux density ratios are invariant 
relative to infrared surface brightness, which suggests that the radiation field that heats 
the cool grains is relatively invariant. However, the 850 fim emission decreases relative to 
70-450 iim emission as the far-infrared surface brightness increases, which suggests that the 
850 /im emission may originate partially from a source other than ~ 25 K dust with a A""^ 
emissivity. 



5. Spectral Energy Distributions 

5.1. Global Spectral Energy Distribution 

To examine the global SED, we will use the global flux density measurements published 
by Dale et al. (2005) with some adjustments. The 3.6 - 8.0 fxm IRAC measurements did not 
include aperture corrections to account for the scattering of light through the detector mate- 
rial; we have applied the "infinite" aperture correction of Reach et al. (2005) to those data. 
The 450 and 850 /xm flux density measurements were repeated for this paper. The SCUBA 
observations cover almost all of the infrared-bright regions in the plane of the galaxy as is 
revealed when comparing the SCUBA images to the 24 /xm images, so a direct measurement 
of the flux densities from the images within the optical disk (a 15.5' x 2.7' eUipse) defined 
in the Third Reference Catalogue of Bright Galaxies (de Vaucouleurs et al. 1991) should be 
accurate. The uncertainties related to missing extended emission will probably be less than 
the the submillimeter calibration uncertainties.^^ 



"^^Note that Dale et al. (2005) applies correction factors of 1.27 and 1.17 to the 450 and 850 /xm observations, 
respectively. The correction factor is calciilated from the fraction of the integrated flux density in the 70 and 
160 /xm bands that falls outside the region covered by the 450 and 850 /xm images. In the case of NGC 4631, 
the 70 and 160 /xm emission falling outside the region observed in the 450 and 850 /im regions is dominated 
by extended emission above and below the plane of the galaxy that may represent the outer regions of the 
PSF from the central structure and other bright sources in the plane of the galaxy. The regions not covered 
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Additionally, we will combine these results with JHK data from the 2MASS Large 
Galaxy Atlas (Jarrett et al. 2003); 12, 60, and 100 /xm integrated flux densities from IRAS 
(Rice et al. 1988); and the 1.23 mm measurements from Dumke et al. (2004) to fill out the 
spectral energy distribution from 1 to 2000 /im. (The 25 fim fiux density from Rice et al. 
(1988), the 870 /im fiux density from Dumke et al. (2004), and the 450 and 850 fim mea- 
surements from Stevens et al. (2005) are discussed in the context of checks on the 24, 450, 
and 850 //m measurements presented here rather than included as redundant flux density 
measurements; see Section 2.) Additionally, CO, free-free, and synchrotron emission (calcu- 
lated in Dumke et al. (2004)) have been subtracted from the 850 and 1230 iim flux densities 
(although the correction is less than 10%). 

The results are presented in Table 3 and Figure 12. This global SED is typical for 
nearby spiral galaxies. Stellar blackbody radiation is seen shortward of 4.5 //m, PAH and 
hot dust emission are seen from 5.7 to 24 /im, and large cool grain emission is seen longward 
of 70 fim. The large cool grain emission has the highest flux density of any source within 
the 1-1000 /im range, as is typical for normal nearby galaxies. 

The slope of the 160-1230 //m emission is shallower than what is expected from black- 
body emission modifled by a X"^ emissivity law for a single temperature component of dust 
(discussed in more detail in Section 5.1.1). This is consistent with the results found by 
Bendo et al. (2003) and Dumke ct al. (2004) although the latter only found excess emission 
at 1230 /im. CO, free-free, and synchrotron emission have already been subtracted from the 
850 and 1230 fira flux densities, so none of these other emission mechanisms can account for 
the observed excess emission. Either a dust component in the 5-15 K range must be present 
or the dust that dominates the 850 - 1230 /xm regime must have an emissivity that does not 
vary as A~^. We examine the possibilities further in the following section. 

As a test of the models, we are searching for models that give gas-to-dust ratios similar 
to the ratio of ~ 165 for the Milky Way, which was calculated by examining the depletion 
of elements in the gas phase of the ISM (Li 2005). Models with signiflcantly higher or 
lower dust mass predictions are probably implausible. This test assumes that the ISM in 
NGC 4631 is similar to the Milky Way's. However, the metaUicity of NGC 4631 is ~ 0.5 
solar metaUicity (Otte et al. 2002). If the gas-to-dust mass ratio varies with metaUicity, then 



in the 450 and 850 /im maps above and below the plane probably do not contain any true sources that 
contributes significantly to the integrated flux density of the galaxy. The 450 and 850 /im PSF is effectively 
truncated in the north-south direction because the observations include chopping to the sky outside the 
plane of the galaxy, so emission from the wings of the PSF would not be present. Therefore, the 450 and 
850 /im flux densities have probably been overcorrected in Dale et al. (2005). We apply no correction to the 
data for the analysis here. 
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the expected ratio for NGC 4631 may be higher than the Milky Way's ratio. 



5.1.1. Simple Modified Blackbody Fitting 

The simplest way to describe the large grain emission longward of 70 /xm is to fit the data 
with one or two blackbodies modified with an emissivity function that scales as a function 
of A"''', where /3 is the emissivity index. This approach is useful just for approximating the 
temperature and emissivity of the bulk of the far-infrared dust emission within the galaxy, 
which can then be used to estimate the mass of the dust. 

The first fit is with a blackbody modified by a emissivity law. This emissivity law 
is what is suggested by simple physical models for far-infrared absorption by small particles 
(e.g. Draine 2004). Figure 13a shows this fit applied to only the 70 - 450 //m data, which 
are the data best fit by this single thermal compoment. The component has a temperature 
of 23 ± 2 K. Excess emission at 850 and 1230 fim is clearly visible in the figure. Note that 
if a single thermal component is fit to the 70 - 1230 fim data, the function significantly 
overestimates the 160 and 450 fim fiux densities while still underestimating the 1230 fim 
flux densities, which is what is expected if the emission at 850 and 1230 fim is in excess 
above the thermal emission expected when extrapolating from shorter wavelengths. 

The mass of this dust can be calculated using 

= 

where D is the distance to the galaxy, f^, is the flux density, Ki, represents the absorption 
opacity of the dust (given by Li & Draine (2001)), and B^{T) is the blackbody function, 
with T the best-fit dust temperature. Using the 450 fim data, we estimate the dust mass 
to be 4.9 X 10^ Mq. For comparison, the total gas mass in this galaxy is 1.4 x 10^° Mq 
(Kennicutt et al. 2003). The gas-to-dust ratio implied by this simple modified blackbody 
fit is 290, which is within a factor of 2 of the Milky Way's ratio. The mass of the hot dust 
is not included in this calculation, but the additional mass from hot dust should bring the 
gas-to-dust ratio closer to the Milky Way's ratio. 

The results show that, to account for the emission at 850 and 1230 fim, either a sec- 
ond thermal component with a cooler temperature is needed or the emissivity law must be 
shallower than A~^. If the source of the excess 850 and 1230 fim emission is a new dust com- 
ponent that produces negligible emission at shorter wavelengths, then this new component 
should not be significantly more massive than the 23 K dust (i.e. the new dust component 
should not increase the dust mass by a factor of 10 or more). The 23 K dust alone can 
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account for at least half of the dust mass expected for this galaxy. 

The next fit is with a single blackbody modified by an emissivity law in which the 
emissivity index is a variable fit to the data. This fit, apphed to the 70 - 1230 iim data, 
is shown in Figure 13b. The fit has a temperature of 28 ± 2 K and an emissivity index of 
1.2 ± 0.1. Note that this accounts for the excess seen at 850 and 1230 //m. However, the 
slope between the 160 and 450 //m data is only marginally consistent with such an emissivity 
function. Since the dust emissivity index is treated as variable, the absorption opacities may 
deviate away from the Li & Draine (2001) values by an uncertain amount. However, if the 
Li & Draine absorption opacities arc still applicable in this fit, the corresponding dust mass 
is 3.3 X 10^ Mq and the gas-to-dust ratio is 420. This is close to the value derived when the 
emissivity index was fixed to -2, although the result with the fixed emissivity are slightly 
closer to the ratio for the Milky Way. 

Another possibility is that the dust should be represented with two blackbodies, each of 
which is modified by emissivity law. To perform this fit properly on dust at temperatures 
near the temperature of the cosmic microwave background (CMB), we cannot use fi, oc 
X~^B^{Tdust) because it may give temperatures colder than Tcmb- Instead, we start with an 
equation that represents the thermal emission from both the CMB and dust in the galaxy 

I, = B,{TcMB)e-^ + B,{Tdust){l - e""), (14) 

where r is the optical depth of the dust. When the background is subtracted, the equation 
is rewritten as 

h{suh) = B,{TcMB)e-^ + B,{Tdust){l - e'^) - B^Tcmb) (15) 

^ r{B,{Tdust) - B,{Tcmb)) t<1 (16) 
Since r is proportional to k,^ (which we are assuming to be proportional to A~^), we can fit 

(X A-2(5,(T,„,,) - B,{Tcmb)) (17) 

to the data. Figure 13c shows the sum of the two thermal components with this modification 
applied to the 70 - 1230 /xm data. The two thermal components have temperatures of 23 ± 2 
and 3.51q8 K. Even though the two components arc fit simultaneously to the 70 - 1230 /xm 
data, the results show that the warmer component represents the majority of the 70 - 450 fim 
emission and that the colder component only contributes a significant amount of emission 
at 850 and 1230 //m. This very cold dust emission is physically implausible because of the 
imphed dust masses. If the dust emissivities in Li & Draine (2001) apply to this cold 
dust component, then the dust mass of the 3.5 K dust is enormous. Based on the fit, 
approximately half of the 850 /im emission comes from this cold dust. Inserting this result 
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into equation 13, we calculate the mass of the 3.5 K dust to be 3.1 x 10^ Mq. This imphes 
a gas-to-dust ratio of ~ 5, which, when compared to the ratio of ~ 165 for the Milky Way, 
is too low to be believable. Therefore, we conclude that 3.5 K dust with an emissivity index 
of -2 cannot give rise to the excess emission observed at 850 - 1230 /im. 

All together, these results show that, using simple blackbodies modified by emissivities 
that scale with wavelength, the excess emission at 850 - 1230 /xm can be represented in two 
ways. The scenario with 3.5 K dust, however, must be rejected because of the implausibly 
high dust masses needed to produce the emission. Further analysis on the SEDs of smaller 
regions within the plane of the galaxy will allow us to further examine whether the variable 
emissivity index scenario is an appropriate descriptions for the 70 - 850 //.m emission. 

5.1.2. Semi- Empirical Model Fitting 

Another option for modeling the dust in this galaxy is to apply the semi-empirical model 
of Dale et al. (2001) to the data. This model is built using a scries of PAH, very small grain, 
and large grain SEDs. These SEDs arc illuminated by a range of radiation fields with relative 
contributions described by the equation 

dMaust oc U-'^dU (18) 

where U is the radiation field. The two variables fit to the data are a and a scahng term. A 
low a indicates relatively strong dust heating, whereas a high a indicates the predominance 
of cooler dust. This is presented as a more complex way to model the dust emission than 
the simple modified blackbody fit. 

The top part of Figure 14 shows the semi-empirical model from Dale et al. (2001) that 
best fits the 5.7 - 850 //m data. This model has an a of 2.31 ± 0.10. This fit does not 
change significantly if the model is fit to only the 3.6 - 850 fim data or the 70 - 850 fim 
data. Of interest here is whether this model describes the dust emission longward of 20 fim. 
The model accurately fits the 70 - 450 /xm data, but it underestimates the 850 - 1230 /xm 
emission. The large grain emission in the Dale et al. (2001) model is represented by a 
series of blackbodies modified with an emissivity function that scales as X~^, where /3 is a 
function of the radiation field (see equation 2 in Dale & Helou 2002). The inabihty of this 
model to reproduce these observed results demonstrates that the summation of a range of 
such modified blackbodies cannot reproduce the excess 850 - 1230 fim emission seen here. 
Either additional dust components or modifications to the emissivity function are needed to 
reproduce the excess submillimeter emission. 

This model can be used to determine the dust mass of the galaxy using the analysis 
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in Section 5.2 of Dale Sz Helou (2002). First, we calculate the dust mass predicted for a 
blackbody modified by a emissivity law with a variable index fit to the 70 - 450 /xm data. 
(We will not use the 850 or 1230 fim. data because the model does not adequately describe 
the emission at those wavelengths.) The best fitting modified blackbody has a temperature 
of 23±2 K and a (3 of 1.8±0.1. Note that this fit statistically differs little from the blackbody 
modified by the emissivity function fit to the 70 - 450 /xm in Section 5.1.1. Next, we use 
the 160 //m flux density and the absorption opacities from Li & Draine (2001) in equation 13 
to calculate a dust mass, which is found to be 5.8 x 10^ Mq. We then apply the appropriate 
correction factor from Figure 6 in Dale & Helou (2002). For a model with a = 2.0 — 2.5, the 
correction factor is ~ 9. This gives the total dust mass as 5.2 x 10^ M©, and the gas-to-dust 
mass ratio as ~ 30. This is still within a factor of five of the ratio of ~ 165 for the Milky Way, 
but it is higher than expected, especially in light of NGC 4631 's metallicity. Nevertheless, 
the dust emission shortward of 450 fim is able to account for all of the expected dust mass. 
This further implies that the modifications that should be applied to this model to explain 
the excess 850 and 1230 /xm emission should not substaintially increase the total dust mass. 



5.1.3. Physical Modeling 

As a third option, we have fit the physical model of Li & Draine (2001) as updated by 
Draine & Li (2006, in preparation) to the global SED. This model is designed to reproduce 
the wavelength-dependent dust extinction and emission within the Milky Way and the Mag- 
ellanic Clouds. The model treats the dust as a mixture of PAHs, amorphous silicate grains 
and carbonaceous grains with a size distribution ranging from tens of atoms to greater than 
1 /im. A fraction 7 of the dust mass is illuminated by a range of radiation flelds described 
by the power law in equation 18 but with the exponent flxed to -2 and the maximum flxed 
to 10^ multiplied by the local interstellar radiation field given by (Mathis et al. 1983). This 
can be thought of as the mass fraction of dust in photodissociation regions. The remainder 
of the dust mass is modeled as a cool cirrus component heated by a single radiation field 
equal to the minimum of the range of radiation fields described by the power-law. Five 
parameters are fit to the data: the minimum value for the illuminating radiation fields, the 
fraction of the total mass in PAH molecules that are smaller than 10^ carbon atoms, the 
starlight continuum intensity in the IRAC bands, 7, and the overall dust mass. This model 
is designed to simultaneously represent the starlight that dominates at 4.5 /im and shorter 
wavelengths, the PAH emission in the the 5.7 and 8 //m wave bands, the 12 and 24 iim hot 
dust emission, and the 60-1230 /im cool cirrus emission. 

The bottom part of Figure 14 shows the physical model that best fits the 1.0-1230 iim 
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data. This physical model is able to reproduce the observed 850 /xm flux densities to within 
the calibration uncertainty of 10%. Furthermore, it is able to account for ~ 75% of the total 
1230 /im emission. In contrast, the single blackbody modified with the emissivity and 
the semi-empirical model both could account for only ~ 30-60% of the 850 fim and 1230 fim 
emission. Unfortunately, the physical model does not accurately reproduce the steep slope 
of the dust emission between 160 and 450 /xm. The model overestimates the 450 /xm flux 
density by ~ 40%, although note that this is only ~ 1.5 times larger than the 25% uncertainty 
in the 450 /im measurement. In summary, it appears that this physical model can flt the 
160-1230 iim data points to within 2.5(7. This could be interpreted as marginally successful, 
although the nature of the mismatch between the model and observations implies that the 
model underpredicts the 1230 /im emission or overpredicts the 450 /im emission. 

Unlike the single thermal component or the semi-empirical models, the physical model 
is able to accurately describe the emission at 850 and 1230 //m because the colors and scale of 
the 160-1230 /im emission can be adjusted independently of the hot dust emission at shorter 
wavelengths. In contrast, the single thermal component model treats the emission between 
70-450 fim as originating from a single temperature when some of the emission at 70 /im 
may be from warmer dust, and the semi-empirical models assume that the colors of all dust 
emission between 5.7 and 1000 fim are interdependent. 

The dust mass of the best-fitting model is 9.7 x 10^ Mq. This suggests a gas-to-dust 
mass ratio of ~ 140, which is very close to the ratio of ~ 165 for the Milky Way. The dust 
mass is only a factor of ~ 2 higher than the dust mass predicted by the 23 K blackbody 
modified by the emissivity given in Section 5.1.1. This is because the temperature of 
most of the dust in this model is ~ 15-18 K, which implies, according to equation 13, a 
mass that is a factor of ~ 2 greater than that calculated using 23 K. Note that this factor 
is significantly lower than the correction factors of ~ 10 — 15 given for the far-infrared or 
submillimeter bands by Figure 6 in Dale & Helou (2002). Also note that the minimum 
temperature predicted by this model is ~ 15 K. A very cold (< 10 K) dust component 
that constitutes approximately 50% or more of the global dust mass may not necessarily be 
needed to explain the 850 or 1230 /xm emission. 

Among the other parameters in the fit, the 7 term (the mass fraction of dust in the 
photodissociation regions) yields the most interesting information. The 7 determined in the 
fit is 0.015 (or 1.5%), which indicates that most of the mass of the dust in the galaxy is 
heated by a weak baseline radiation field in diffuse cirrus emission. The model also predicts 
that the small fraction of dust in photodissociation regions produces ~ 10% of the total dust 
emission. 

The other parameters from the physical model fit are the PAH mass fraction, which is 
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0.034, and the minimum value of the illuminating radiation field, which is 2.0 times the local 
interstellar radiation field. 



5.2. Spectral Energy Distributions in Discrete Regions 

The SEDs of discrete regions within this galaxy will reveal where the excess 850 //m 
emission is present in the galaxy. This will provide additional clues on where the submil- 
limeter excess arises from, which can then be used to better understand the nature of the 
excess. We composed SEDs for individual regions throughout the disk using the data for 
the regions in Table 2 but excluding locations where either 450 or 850 /xm data were not 
available. Below, we will attempt to describe the data with both simple and complex models 
of the emission. 



5.2.1. Simple Modified Blackbody Fitting 

First, we will fit the 70-450 data with a single blackbody modified by a emis- 
sivity law. Example plots of the fits for the central structure, bright knots in the disk, and 
faint regions in the disk are shown in Figure 15. The parameters for these fits and the 
difference between the measured 850 fim flux density and the fits (henceforth referred to as 
the excess 850 //.m emission) is presented in Table 4. Additionally, Figures 16 and 17 show 
how the temperature and excess 850 /im emission from these fits vary with the total infrared 
luminosity (calculated using the MIPS 24, 70 and 160 //m data and equation 4 in Dale & 
Helou 2002). 

These results show that the function fit to the 70 - 450 /xm data represents those data 
very well. As seen in Figure 16, very little variation in the temperature is present, as 
was expected from the invariance of the 70/160 /xm and 160/450 //m colors in Section 4.4. 
Additionally, these results provide us with two pieces of evidence that suggest that the excess 
850 yum emission cannot be explained simply by an emissivity function that deviates from 
the law given by Li & Draine (2001). First are the results from the fits to the 70 - 
450 fim data that show that dust with an emissivity index of -2 is adequate for describing 
the dust emission. Second is the elbow-like bend in the SEDs from 160 to 850 /xm found 
in many of the SEDs. This bend cannot be produced by varying the emissivity index for 
the 70-1000 iim dust emission. If the emissivity index is variable, it only varies longward of 
450 /im. It is possible, however, that some constituent of the interstellar medium other than 
large solid grains contributes to the emission at 850 /im. 
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The variations in the excess 850 //m emission in Figure 17 show that the excess becomes 
more prominent in the 850 /xm emission from infrared-faint regions. In the infrared-bright 
regime, the contribution of the excess to the total emission is only ~ 20%. This is close to 
the calibration uncertainties from SCUBA as well as other similar submillimeter data. The 
excess therefore becomes statistically indistinguishable from the ~ 25 K dust emission found 
within the waveband. In infrared-faint regions, however, the excess 850 /xm emission can 
contribute ~ 80% of the total 850 fj,m flux density. 

Note that the excess 850 /xm emission observed here (and in Sections 5.2.2 and 5.2.3) as 
well as the variations in the 70/850, 160/850, and 450/850 /xm flux density ratios is unlikely 
to be produced by either a background offset or a calibration error. If the background were 
incompletely subtracted, the excess 850 /tm emission listed in Table 4 would be constant. 
However, the excess 850 /im emission instead varies signiflcantly, with low excess values 
generally more likely to be found in infrared- faint regions (although the trend between the 
850 /im excess emission and the infrared luminosity is statistically weak). If the calibration of 
the 850 /xm emission was incorrect, then the trend in Figure 17 would be flat. Moreover, since 
the 450 and 850 //.m data were processed using the same software, calibration standards, and 
background regions, any offset in the 850 fim data should also manifest itself in the 450 fim 
data. Since the 450 /im data are consistent with expectations based on shorter wavelengths 
but the 850 /tm data are not, we conclude that the excess 850 /tm emission is real and not 
an artifact of the data processing. 



5.2.2. Semi- Empirical Model Fitting 

Next, we have flt the 5.7 - 450 /xm data with the semi-empirical model of Dale et al. 
(2001). Example plots of these flts are shown in Figure 18, and parameters from the flts are 
shown in Table 5. We did attempt fits that included the 850 /tm data points, but the match 
to the 850 /tm data did not improve significantly despite its inclusion, and the fits failed to 
converge on appropriate a values in some situations. We therefore excluded the 850 /tm data 
from the fits. 

As with the single blackbodies modified with the emissivity laws, these fits still 
reveal the presence of excess 850 /tm emission. Plots of how the excess 850 /tm emission 
derived from these model fits vary with total infrared luminosity are shown in Figure 19. 
Again, the fraction of 850 /tm flux density from the excess increases as the total infrared 
luminosity decreases. In terms of determining the presence and magnitude of the excess 
850 /tm emission, few differences exist between using these semi-empirical model and using 
single blackbodies modified with emissivity laws. 
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The a parameter does show a correlation with total infrared luminosity, as shown in 
Figure 20. Given that the 70/160 /xm flux density ratio does not vary for these regions but 
the 24/160 fim docs vary, the increase in heating must be reflected primarily in changes in 
the amount of hot dust emission. Again, this is somewhat consistent with the results of 
Dale et al. (2001), which showed that the 24 fim emission would vary most strongly with 
dust heating. The large grain emission responsible for the 70 and 160 /xm emission, however, 
must be relatively unaffected. 

5.2.3. Physical Modeling 

Finally, we present the physical models that have been fit to the 3.6 - 850 /im data. 
Example plots of these fits arc shown in Figure 21. The parameters from the fits arc given in 
Table 6 and the 450 and 850 /xm excess flux densities are given in Table 7. Unlike the single 
blackbody fits and the semi-empirical model fits, these physical models perform much better 
at describing the 850 /im data. However, a significant fraction of the 850 /im fiux density 
(up to 30%) can still be identified as in excess above the model predictions. Moreover, 
as seen in Figure 22, this ratio of the excess to total 850 /xm fiux density still varies with 
the infrared luminosity. Table 7 also shows that the model overestimates the 450 /xm fiux 
density in the same infrared- faint regions where it underestimates the 850 ^m emission. (The 
overestimates are expressed as negative 450 fim excess emission.) The overestimates of the 
450 //.m emission may reach up to ~ 80% in some locations. Nonetheless, these mismatches 
only reach a maximum of ~ 3cr significance within individual regions. So, we conclude that 
these physical models perform marginally well in explaining the 450 and 850 /im emission 
but that excess 850 //m emission may still be present. 

The PAH mass fraction determined for the 40" subregions do show a variation with 
radius. This is shown in Figure 23. The monotonic decrease in the PAH mass fraction with 
radius is similar to the variation with radius seen in the (PAH 8)/160 /im fiux density ratio 
in Figure 9. Again, no local peaks in the PAH mass fraction corresponding to local peaks in 
the total infrared surface brightness are seen outside the center of the galaxy, at least on the 
kiloparsec scales of the 40" regions. This additional analysis demonstrates that variations in 
PAH emission in this galaxy are a function of distance from the center, not infrared surface 
brightness. 

Also of interest are the variations in 7 (the mass fraction of dust in photodissociation 
regions) with radius. This is shown in Figure 24. Note that this parameter stays at the level 
of 0.01 - 0.02 except in the central structure, where it reaches up to ~ 0.03. This shows 
that the relative mass of dust in photodissociation regions is twice as high as in the central 



-27- 



structure, although it is still a relatively small amount of the total dust mass. 

5.3. Discussion 

In summary, the SED results show that the far-infrared / submillimeter dust emission 
is somewhat unusual in two respects. First is the invariance in the color temperature of dust 
emission in the 70 - 450 /im range. This result is unusual in comparison to other recent 
Spitzer observations of galaxies. However, this invariance could be explained if the bulk of 

the dust emission originates from dust outside of regions with strong heating by young stars, 
which may be consistent with the relatively low mass fraction of dust in photodissociation 
regions found when the physical model was fit to the global SED. The increase in infrared 
luminosity would then primarily be tracing an increase in dust mass integrated along the line 
of sight rather than an increase in dust heating. Note that, because this is an edge-on galaxy, 
increases in the integrated dust mass can be achieved if the line of sight passes through either 
a greater fraction of diffuse dust distributed evenly throughout the galaxy's disk or massive 
individual dust clouds within the galaxy. 

The other unusual result from these observations is the slope of the dust emission be- 
tween 160 and 1230 //m. The single modified blackbody fits and the semi-empirical model 
fits reveal the presence of 850 and 1230 /im emission in excess of what is anticipated from 
models of thermal dust emission. The physical models describe the data better, but the 
results from the fits imply that such an excess at 850 and 1230 /im may be present. 

The characteristics of the source of the excess emission - if it is present - are complex. 
First, if it originates from dust, then the modifications that need to be applied to the models 
to explain the excess emission should not substantially increase the total estimated dust 
mass. This is because both simple modified blackbodies and more sophisticated models fit 
to the 5.7 - 450 fim emission give dust masses that suggest dust-to-gas mass ratios similar 
to what is measured in the Milky Way. Second, if the excess emission represents variability 
in the emissivity index of the dust emission, the variation only occurs longward of 450 /xm, 
as explained in Section 5.2.1. Third, the 850 //.m excess appears to be found throughout the 
disk. However, if the excess is present in the brightest infrared regions, it is relatively weak 
compared to 20-25 K large grain dust emission at 850 //m. Finally, the 850 //m excess also 
diminishes in strength in infrared-faint regions, although it may still contribute 30-80% of 
the 850 /im emission. 

Klaas et al. (2001), Dunne & Eales (2001), and Bendo et al. (2003) had reported a 
degeneracy in trying to represent the far-infrared to submillimeter dust emission either as 
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one blackbody modified by an emissivity power law with a variable index or as multiple 
blackbodies modified with fixed X"^ emissivity laws. The results here have shown that 
the variable emissivity index scenario cannot accurately describe the 70-850 /im emission of 
discrete regions within a galaxy, so the approach should not be applied globally. The 850 ^m. 
excess is best described as a component of the SED that is separate from the thermal dust 
emission seen at 70-450 fim.. However, the 850 /xm excess should not be described as a very 
cold thermal dust component because the implied dust masses are too high. 

Some previous observations, such as Regan ct al. (2004), may have missed similar excess 
850 ^m. emission from nearby galaxies because of the nature of the observations. Most global 
SEDs will be dominated by a few bright regions. If SEDs for subregions are extracted, the 
subregions are usually of the brightest sources within the targets. The results for this galaxy 
demonstrate that excess submillimeter and millimeter emission is small compared to the 
~ 20 K large grain emission and cannot be easily detected. However, the excess becomes 
easier to detect in fainter regions in the disk, as it may constitute 80% of the emission in 
some regions. Therefore, we suggest that the search for excess submillimeter or millimeter 
emission in other galaxies should not be performed with global or nuclear SEDs but instead 
should be carried out in fainter regions with galaxies' optical disks. 

Through a process of elimination, only a couple of options are available for describing 
the 850 and 1230 /xm excess if it is indeed present. One possibility is that the excess emission 
originates from dust grains with high submillimeter emissivities relative to their absorption 
cross sections. Such grains, such as the fractal grains proposed by Reach et al. (1995) for the 
Milky Way and Dumke et al. (2004) for NGC 4631, may be able to maintain relatively low 
temperatures, so they would be only detectable at submillimeter wavelengths. This scenario, 
however, may be physically implausible (Li 2005). Another possibility, also considered by 
Reach et al. (1995) for the Milky Way, is the possibility that the emissivity of the ~ 20-25 K 
large grains is enhanced above what is expected by Li & Draine (2001) by resonances in the 
dust grains caused by impurities. Further physical modeling as well as more observations of 
submillimeter excesses will be needed to better understand this phenomenon. 

6. Conclusions 

The primary results from this research can be summarized as follows: 

1. PAH emission at 8 /xm is closely related to hot dust emission at 24 /xm. The relation 
holds on 1.7 kpc scales but begins to break down at 650 pc scales. Furthermore, the PAH 
emission is not as centrally peaked as the hot dust emission on scales of hundreds of pc. 
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2. Variations in the strength of PAH emission relative to cool dust emission and in the 
fraction of dust mass in PAHs depend primarily on radius, not infrared surface brightness. 

3. The 70 - 450 //m color temperature does not appear to vary with surface brightness in 
this galaxy, even though the 24/70 //m and 24/160 //m flux density ratios do vary with surface 
brightness. This implies that a substantial part of the dust emission in this wavelength regime 
may originate from a cool, diffuse cirrus component. 

4. The 850 and 1230 //m emission in this galaxy is found to exceed what is anticipated 
from either the 23 K blackbody modified by a emissivity law or the semi-empirical model 
that describes the emission shortward of 850 /xm. The physical dust models marginally 
describe the data, but they leave open the possibility that excess emission at 850 /im is 
present. If present, the 850 excess is highest in regions of moderate infrared brightness 
and constitutes the greatest fraction of the 850 jim emission in infrared-faint regions. 

Two major implications for using PAH emission as a star formation indicator arise 
from this research. First, these results demonstrate that the spatial correlation between 
PAH emission and 24 fim dust emission breaks down on scales of hundreds of parsecs. This 
implies that PAH emission cannot be used as an accurate star formation tracer on such scales 
(if 24 iim emission primarily traces star formation regions). Second, these results suggest 
that the ratio of PAH to dust mass varies radially, although the reason for the radial variation 
needs further study. If PAH emission is to be used as a tracer of dust or star formation, 
then these radial effects need to be taken into account. Further comparisons using 8, 24, and 
160 fjim data of nearby galaxies are needed to understand how PAH emission relates to dust 
emission in other galaxies and to understand the radial variations in the PAH emission. 

The color variations and SEDs found within NGC 4631 have some major implications 
for dust modeling. The absence of variations in the 70/160 //m and 160/450 iim flux density 
ratios in this galaxy but the presence of variations in the 70/160 /xm ratio in other galax- 
iesneeds to be studied further so as to understand the conditions that cause the ratios to 
vary or remain constant. Second, dust models may need to be adjusted to better account for 
the excess emission at 850 fim if it is present. Further analysis with Spitzer 70 and 160 fim 
data as well as 850 /im data with comparable resolution are needed not only to confirm the 
presence of the 850 /im excess emission in other galaxies but to provide better descriptions 
of where the excess can be found. Furthermore, additional observations of nearby galaxies 
at multiple wavelengths between 160 and 850 /im are needed so as to better constrain the 
shape of the SED in this wavelength regime. 
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Fig. 1. — Global images of NGC 4631 from 3.6 to 850 fim. Each image is 18' x 6', with north 
up and east to the left. Circles representing the FWHM of the PSF have been plotted in the 
lower left corner of the 24 - 850 fim images. In the 8 /im image, a muxbleed effect is visible 
as streaking to the northeast and southwest of the center of the galaxy. Also note that if 
the stellar continuum is subtracted from the 8 /im image using equation 1, the resulting 
PAH emission map is superficially similar to the above 8 /^m image. In the 70 /im image, 
a latent image effect has created negative streaking to the south of the center. This latent 
image effect also results in an artifact in the background subtraction that creates the positive 
streaking to the north of the center. The edges of the usable parts of the 450 and 850 /im 
images is outlined in green. For display purposes only, the 450 /im image has been smoothed 
with an 8" Gaussian, and the 850 /xm image has been smoothed with a 15" Gaussian. 
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Fig. 1. — Continued. 
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Fig. 2. — Images of the center of NGC 4631 from 3.6 to 850 /xm. Each image is 3' x 1'.5, with 
north up and east to the left. For display purposes only, the 450 fim image has been smoothed 
with an 8" Gaussian, and the 850 fim image has been smoothed with a 15" Gaussian. The 
edge of the usable part of the 450 mum observations is outlined in green. 
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Fig. 3. — Global images of NGC 4631 at 24 jjm with the 15" and 40" extraction regions 
overlaid as circles on the image. The cyan circles mark the regions plotted in Figures 15, 18, 
and 21. Each image is 18' x 6', with north up and east to the left. 
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Fig. 4. — A comparison of the PAH 8 /um and 24 jj,m flux densities measured in the 15" 
(650 pc) regions listed in Table 1 and the 40" (1.7 kpc) regions listed in Table 2. Uncertainties 
from background noise are smaller than the symbols. The lines show the best fitting relations 
between the (PAH 8)/24 //m flux density ratio and the 24 //m flux density. 
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Fig. 5. — Global image of the (PAH 8)/24 fim flux density ratio in NGC 4631. This image is 
14' X 3', with north up and east to the left. The PAH 8 /xm data were convolved with a kernel 
described in Section 2.4 to match the resolution of the 24 /xm data. Blue colors correspond 
to regions with relatively strong PAH 8 fim emission; red colors correspond to regions with 
relatively strong 24 /im emission. Pixels with low signal-to-noise ratios have been left blank, 
and the regions to the northeast and southwest of the center that are heavily affected by 
muxbleed in the 8 /im image have been masked out. The color bar on the right shows the 
correspondence between the shading and the numerical values for the (PAH 8)/24 /im flux 
density ratio. 
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Fig. 6. — PAH 8 /iin surface brightness radial profiles (shown as green dashed lines with 
uncertainties shown in light green) and 24 fim surface brightness radial profiles (shown as 
red dotted lines with uncertainties shown in pink) of ten regions in the disk of NGC 4631. 
All PAH 8 fim data are convolved with the kernels described in Section 2.4 before radial 
profiles are extracted. The radial profiles are all normalized to 1 at the innermost radii. 
The region centers are identified by J2000 right ascension and declination coordinates in the 
lower left corner of each plot. Regions where uncertainties overlap are shown in dark yellow. 
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Fig. 7. — A comparison of the 24 and 70 /xm flux densities measured in the 15" regions hsted 
in Table 1, and a comparison of the 24 and 160 fim flux densities measured in the 40" regions 
listed in Table 2. Uncertainties from background noise are smaller than the symbols. The 
lines show the best fitting relations between the flux density ratios on the y-axes and the 
flux densities on the x-axes. The two data points that deviate from the relation in the top 
plot are marked with circles. These data points are discussed in the text. 
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Fig. 8. — A comparison of the PAH 8 and 160 /im flux densities measured in the 40" regions 
hsted in Table 2, and a comparison of the PAH 8 and 850 //m flux densities measured in 
the 15" regions listed in Table 1. Uncertainties from background noise are smaller than the 
symbols except where error pars are shown. The lines show the best fitting relations between 
the flux density ratios on the y-axes and the flux densities on the x-axes. 
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Fig. 9. — A plot of the total infrared luminosity (top), the 24/160 /xm ratio (middle), and 
the (PAH 8)/160 //m ratio (bottom) plotted as a function of distance from the center. These 
data are calculated using the flux densities measured in the 40" regions listed in Table 2. 
The two regions listed in Table 2 that are located off the major axis of optical disk of the 
galaxy are excluded from this plot. Uncertainties from background noise in the middle and 
bottom plots are smaller than the symbols. Data from the east side of the galaxy are plotted 
to the left, and data from the west side are plotted to the right. 
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Fig. 10. — A comparison of the 70, 160, and 450 /im flux densities measured in the 40" regions 
listed in Table 2. Except where uncertainty bars are shown, uncertainties from background 
noise are smaller than the symbols. The hnes show the best fitting relations between the 
flux density ratios on the y-axes and the flux densities on the x-axes. 
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Fig. 11. — A comparison of the 70 and 850 //m flux densities measured in the 15" regions 
hsted in Table 1, and comparisons of the 160, 450 //m, and 850 /xm flux densities measured 
in the 40" regions Ustcd in Table 1. Except where uncertainty bars are shown, uncertainties 
from background noise are smaller than the symbols. The lines show the best fitting relations 
between the flux density ratios on the y-axes and the flux densities on the x-axes. 
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Fig. 12. — The 1.3 - 1230 jj^m global SED for NGC 4631. Except where uncertainty bars 
are shown, the uncertainties are smaller than the symbols in this plot. Note that the 25 /xm 
measurement from Rice et al. (1988), the 870 /xm measurement from Dumke et al. (2004), 
and the 450 and 850 /xm measurements from Stevens et al. (2005) are not included in this 
plot, although they are discussed in the text. 
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Fig. 13. — The 30 - 1230 /xm global SED for NGC 4631 with various simple modified 
blackbodies fit to the data. Except where the uncertainty bars are shown, the uncertainties 
are smaller than the symbols in this plot. The functions are fit to the 70 - 1230 iim data 
except where noted. Part a features one blackbody modified with a emissivity law (fit 
only to the 70 - 450 fim data), with the best fit having a temperature of 23 ± 2 K. Part b 
features one blackbody modified with a emissivity law with an index that was fit to the 
data. The best fit has a temperature of 28 ± 2 K and a of 1.2 ± 0.1. Part c features two 
blackbodies modified with emissivity laws (with an adjustment to account for the CMB; 
see the text for details). The individual thermal components are shown as dotted lines. The 
best fitting functions had temperatures of 23 ±2 and 3.5^o!8 
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Fig. 14.— The 1.3 - 1230 //m global SED for NGC 4631 with the best fitting semi-empirical 
dust emission model of Dale et al. (2001) (with a of 2.32 ± 0.10) overlaid on the data in 
panel a and the best fitting physical dust emission model overlaid on the data in panel b. 
See the text for the parameters for the model in part b. Except where uncertainty bars are 
shown, the uncertainties are smaller than the symbols in this plot. 
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Fig. 15. — 70 - 850 /im SEDs for example 40" regions listed in Table 2 with the best fitting 
blackbody modified by a emissivity law overlaid on the data. The top row consists of 
data from regions in the central structure. The middle row consists of data from three bright 
knots in the outer disk. The bottom row consists of data from three faint regions in the outer 
disk. The J2000 coordinates of the regions are listed in the top right corners of each plot. 
Except where uncertainty bars are shown, the uncertainties are smaller than the symbols in 
this plot. 
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Fig. 16. — A comparison of the temperatures of blackbodies modified with emissivity 
laws fit to the 70 - 450 /xm data to the total infrared luminosities for the 40" regions listed 
in Table 2. 



-50- 




Total Infrared Luminosity {L^ 



Fig. 17. — A comparison of the ratio of the excess 850 fim flux density to total 850 fim flux 
density to the total infrared luminosities for the 40" regions listed in Table 2. The excess 
emission in this plot is derived from the difference between the total 850 /im emission and 
the 850 /im emission expected from a blackbody modified with a emissivity law fit to 
the 70 - 450 //m data. 
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Fig. 18. — 3.6 - 850 /im SEDs for example 40" regions listed in Table 2 with the semi- 
empirical models of Dale et al. (2001) that best fit the 5.7 - 450 /xm overlaid on the data. 
These arc the same regions shown in Figure 15. The J2000 coordinates of the regions are 
listed in the top left corners of each plot. Except where uncertainty bars are shown, the 
uncertainties are smaller than the symbols in this plot. 
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Fig. 19. — A comparison of the ratio of the excess 850 fim flux density to total 850 fim flux 
density to the total infrared luminosities for the 40" regions listed in Table 2. The excess 
emission in this plot is derived from the difference between the total 850 /im emission and 
the 850 iim emission expected from the semi-empirical models of Dale et al. (2001) fit to the 
5.7 - 450 iim data. 
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Fig. 20. — A comparison of the a of the semi-empirical models fit to the 5.7 - 450 /xm data 
to the total infrared luminosities for the 40" regions listed in Table 2. Lower a represent 
regions with proportionally more heating by stronger radiation fields. 
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Fig. 21. — The same as Figure 18 but with the best fitting physical models overlaid on the 
data. 
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Fig. 22. — A comparison of the ratio of the excess 850 fim flux density to total 850 fim flux 
density to the total infrared luminosities for the 40" regions listed in Table 2. The excess 
emission in this plot is derived from the difference between the total 850 /im emission and 
the 850 pLUi emission expected from the physical models fit to the 3.6-850 iim data. The 
dashed line demarcates where the y-axis equals 0. 
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Fig. 23. — The PAH mass fraction determined for the 40" regions hsted in Table 2 plotted 
as a function of distance from the center. Data from the east side of the galaxy are plotted 
to the left, and data from the west side are plotted to the right. 
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Fig. 24. — The 7 term determined for the 40" regions Usted in Table 2 plotted as a function 
of distance from the center. This 7 term represents the fraction of the dust mass illuminated 
by the radiation fields described by the power law (i.e.the dust heated in photodissociation 
regions). Data from the east side of the galaxy are plotted to the left, and data from the 
west side are plotted to the right. 
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.0 


+32 


32 


50 


2.392 X 


10" 




1.575 X 10" 




5.167 X 


10" 




1.3298 X 10" 


-2 


1.2519 X 10" 




1.667 X 10-1 


1.00 X 


10" 




12 42 


20, 


.7 


+32 


31 


46 


8.41 X 10~ 


4 


5.66 X 10- 


4 


1.173 X 


10" 


-3 


2.798 X 10- 


3 


3.791 X 10- 


3 


7.13 X 10-^ 


4.6 X 


10" 


-3 


12 42 


21 


.4 


+32 


33 


05 


1.645 X 


10" 


-3 


1.157 X 10" 


-3 


4.135 X 


10" 


-3 


1.1073 X 10" 


-2 


1.5873 X 10" 


-2 


1.659 X 10-1 


8.6 X 


10" 


-3 


12 42 


21 


.5 


+32 


32 


23 


1.882 X 


10" 


-3 


1.246 X 10" 


-3 


3.499 X 


10" 


-3 


8.500 X 10" 


3 


9.754 X 10- 


3 


1.467 X 10"! 


8.2 X 


10" 


'3 


12 42 


22, 


.0 


+32 


32 


44 


2.926 X 


10" 


-3 


2.067 X 10" 


-3 


7.723 X 


10" 


-3 


2.0298 X 10" 


-2 


3.3467 X 10" 


-2 


3.708 X 10"! 


1.65 X 


10 


-2 


12 42 


23, 


.2 


+32 


31 


44 


5.83 X 


10" 


4 


4.07 X 10" 


4 


7.78 X 


10" 


4 


2.070 X 10" 


3 


2.687 X 10- 


3 


5.82 X 10"^ 


5.6 X 


10" 


-3 


12 42 


23, 


.5 


+32 


32 


44 


2.036 X 


10" 


-3 


1.417 X 10" 


-3 


4.648 X 


10" 


-3 


1.1665 X 10" 


-2 


1.4880 X 10" 


-2 


2.060 X 10-1 


1.35 X 


10" 


-2 


12 42 


23 


.6 


+32 


32 


25 


1.491 X 


10" 


-3 


9.89 X 10- 


4 


2.432 X 


10" 


-3 


6.961 X 10- 


3 


7.083 X 10- 


3 


1.189 X 10-1 


1.19 X 


10" 


-2 


12 42 


24, 


.0 


+32 


32 


08 


1.133 X 


10" 


-3 


7.67 X 10- 


4 


1.676 X 


10" 


-3 


3.909 X 10- 


3 


5.996 X 10- 


3 


1.187 X 10-1 


7.4 X 


10" 


-3 


12 42 


24, 


.8 


+32 


32 


14 


1.211 X 


10" 


-3 


8.21 X 10- 


4 


1.717 X 


10" 


-3 


4.296 X 10- 


3 


6.017 X 10- 


3 


1.166 X 10-1 


1.07 X 


10" 


-2 


12 42 


25, 


.3 


+32 


32 


50 


1.377 X 


10" 


-3 


9.21 X 10- 


4 


2.344 X 


10" 


-3 


6.977 X 10- 


3 


5.666 X 10- 


3 


9.11 X 10-^ 


1.12 X 


10" 


-2 


12 42 


26, 


.0 


+32 


32 


31 


1.098 X 


10" 


-3 


7.39 X 10- 


4 


1.423 X 


10" 


-3 


3.563 X 10- 


3 


3.900 X 10- 


3 


7.69 X 10-^ 


6.1 X 


10" 


-3 


12 42 


27, 


.3 


+32 


32 


21 


8.46 X 10" 


4 


5.64 X 10- 


4 


9.54 X 


10- 


4 


2.553 X 10- 


3 


2.687 X 10- 


3 


6.94 X 10-^ 


2.7 X 


10" 


-3 


12 42 


27, 


.6 


+32 


32 


41 


9.95 X 


10- 


4 


6.62 X 10- 


4 


1.188 X 


10" 


-3 


3.017 X 10- 


3 


3.343 X 10- 


3 


6.72 X 10-^ 


6.6 X 


10" 


-3 


12 42 


28, 


.7 


+32 


32 


40 


9.94 X 


10- 


4 


6.66 X 10- 


4 


1.246 X 


10" 


-3 


3.322 X 10- 


3 


3.894 X 10- 


3 


7.37 X 10-^ 


1.38 X 


10" 


-2 



00 



Table 1 — Continued 



Position (J2000) 








Flux Density (Jy)'* 












R.A. 


Dec. 


3.6 iJ.in^ 


4.5 /Lim^ 


5.7 iJ.m'^ 


8.0 /am« 


24 fim' 


70 (imS 




850 




12 42 29.3 


+32 32 49 


8.25 X 10~* 


5.59 X 10-* 


1.043 X 10-^ 


2.833 X 10-^ 


3.644 X 10-^ 


6.66 X 10- 


2 


8.6 X 


10-^ 


12 42 30.8 


+32 32 31 


6.98 X lO-** 


4.94 X lO-** 


7.80 X 10-* 


2.277 X 10-^ 


4.414 X 10-^ 


8.35 X 10- 


2 


4.8 X 


10-^ 


12 42 32.6 


+32 32 37 


4.48 X 10-" 


3.25 X 10-" 


2.61 X 10-" 


1.105 X 10-^ 


1.485 X 10-^ 


3.34 X 10- 


2 







^The 3.6 - 24 /Ltm flux densities are measured in maps that have been convolved with the kernels described in Section 2.4 so as to match 
the higher-resolution PSFs to the PSF of the MIPS 70 fxm data. 

^Calibration uncertainty: 10%; background noise uncertainty: 4 X 10~^ Jy 

"^Calibration uncertainty: 10%; background noise uncertainty: 6 X 10~^ Jy 

•^Calibration uncertainty: 10%; background noise uncertainty: 1.7 X 10~^ Jy 

^Calibration uncertainty: 10%; background noise uncertainty: 1.4 X 10~^ Jy; note that these data are not stellar continuum subtracted. 
^Calibration uncertainty: 10%; background noise uncertainty: 1.1 X 10~^ Jy 
^Calibration uncertainty: 20%; background noise uncertainty: 7 X 10~^ Jy 

^Calibration uncertainty: 10%; background noise uncertainty: 1.0 X 10~^ Jy; note that CO, free-free, and synchrotron emission have not 
been subtracted from these data, but since the sum of these emission components only contributes less than 10% to the global flux density, 
the absence of a correction for these data should not strongly affect the results. 



Table 2. Flux Densities Measured within Discrete 40" Regions in NGC 4631 





Pc 

R 


V, 


ion 


r,T2()0()) 




:!.(, ,,m'' 




4,."i / 






- - -1 
,1111 




/(111' 


Flii 


24 ,,iii' 


■v 


70 / 




KiO /mi'' 


4.".0 






,1111 


i 


12 


41 


37, 


6 


+ 32 


:ii 


34 


1,94 X 10^ 


3 


1,00 X 


10" 


3 


2.48 X 10" 


3 


5.49 X 10" 


3 


9.8,'') X 10" 


3 


1.849 X 


10"! 


4,72 X 10"! 






3.7 X 


10" 


2 


12 


41 


40. 


,9 


+ 32 


31 


50 


3,95 X 10" 


3 


2.28 X 


10" 


3 


5.87 X 10" 


3 


1.349 X 10" 


-2 


2.945 X 10" 


2 


4.938 X 


10"! 


1,236 X 10° 


1.66 X 


10-1 


4.7 X 


10- 


2 


12 


41 


43. 


,2 


+32 


31 


56 


4.64 X 10~ 


3 


2.60 X 


IQ- 


3 


7.25 X 10- 


3 


1.667 X 10" 


-2 


3.064 X IQ- 


2 


5.232 X 


10-1 


1.383 X 10° 


1.65 X 


10-1 


6.4 X 


10- 


2 


12 


41 


45. 


,8 


+32 


32 


00 


6.29 X 10~ 


3 


3.77 X 


10- 


3 


1.177 X 10" 


-2 


2.739 X 10" 


-2 


4.964 X 10- 


2 


7.647 X 


10-1 


2.070 X 10° 


1.77 X 


10-1 


7.1 X 


10- 


2 


12 


41 


49. 


,4 


+32 


32 


04 


9.49 X 10~ 


3 


6.11 X 


10- 


3 


1.969 X 10" 


-2 


4.593 X 10" 


-2 


8.639 X 10- 


2 


1.1708 


X 10° 


3.294 X 10° 


2.48 X 


10-1 


9.1 X 


10- 


2 


12 


41 


51. 


,8 


+32 


32 


10 


1.083 X 10" 


-2 


6.74 X 


10- 


3 


2.211 X 10" 


-2 


6.146 X 10" 


-2 


8.308 X 10- 


2 


1.2284 


X 10° 


3.242 X 10° 


2.72 X 


10-1 


9.5 X 


10- 


2 


12 


41 


54, 


2 


+32 


32 


14 


1,497 X 10" 




9.35 X 


10- 




2.841 X 10" 




6.537 X 10" 




1.1305 X 10" 




1.5426 


X lO" 


3.967 X 10° 


4.22 X 


10-1 


1.05 X 


10" 




12 


41 


58, 


,4 


+32 


32 


12 


2,153 X 10" 


-2 


1.328 X 


10" 


-2 


3.617 X 10" 


-2 


8.042 X 10" 


-2 


9.937 X 10" 


2 


1.4413 


X lo" 


4.006 X 10° 


5.19 X 


10"! 


9.5 X 


10" 


2 


12 


42 


01, 


,3 


+32 


32 


17 


3,402 X 10" 


-2 


2.141 X 


10" 


-2 


6.719 X 10" 


-2 


1.5801 X 10 


-1 


2.0379 X 10" 


-1 


2.3387 


X lo" 


6.470 X 10° 


7.15 X 


10"! 


1.21 X 


10" 


-1 


12 


42 


04, 


,3 


+32 


32 


25 


7,091 X 10" 


-2 


4.637 X 


10" 


-2 


1.5671 X 10 


-1 


3.8625 X 10 


-1 


5.6175 X 10" 


-1 


5.4659 


X 10° 


1.4361 X K)! 


1.332 


X lo" 


2.11 X 


10" 


-1 


12 


42 


07, 


,9 


+32 


32 


35 


9.913 X 10" 


-2 


6.584 X 


10" 


-2 


1.9151 X 10 


-1 


4.7155 X 10' 


-1 


7.2293 X 10" 


-1 


6.2692 


X 10° 


1.6281 X lOl 


1.457 X lO"^ 


2.16 X 


10" 


-1 


12 


42 


10, 


,9 


+32 


32 


35 


6.881 X 10" 


-2 


4.729 X 


10" 


-2 


1.5475 X 10 


-1 


3.8830 X 10 


-1 


6.6790 X 10" 


-1 


5.4380 


X 10° 


1.5191 X IqI 


1.375 


X 10° 


1.98 X 


10" 


-1 


12 


42 


14, 


,0 


+32 


32 


38 


3.300 X 10" 


-2 


2.198 X 


10" 


-2 


6.623 X 10" 


-2 


1.6114 X 10 


-1 


1.9198 X 10" 


-1 


2.4416 


X 10° 


6.474 X 10° 


8.26 X 


10-1 


1.10 X 


10" 


-1 


12 


42 


16, 


,7 


+32 


32 


39 


1.850 X 10" 


-2 


1.200 X 


10" 


-2 


3.206 X 10" 


-2 


7.759 X 10" 


-2 


8.832 X 10- 


2 


1.2228 


X 10° 


3.664 X 10° 


6.55 X 


10-1 


7.6 X 


10- 


2 


12 


42 


19, 


,7 


+32 


32 


37 


1.409 X 10" 


-2 


9.20 X 


10- 


3 


2.361 X 10" 


-2 


6.692 X 10" 


-2 


7.638 X 10- 


2 


1.0366 


X 10° 


3.002 X 10° 


3.98 X 


10-1 


7.1 X 


10- 


2 


12 


42 


20, 


,2 


+32 


31 


59 


7.91 X 10- 


3 


6.23 X 


10- 


3 


1.071 X 10" 


-2 


2.518 X 10" 


-2 


3.669 X 10- 


2 


6.263 X 


10-1 


1.632 X 10° 


1.42 X 


10-1 


3.2 X 


10- 


2 


12 


42 


22, 


,0 


+32 


32 


44 


1.288 X 10" 


-2 


8.71 X 


10- 


3 


2.323 X 10" 


-2 


6.746 X 10" 


-2 


9.974 X 10- 


2 


1.2033 


X 10° 


3.196 X 10° 


3.73 X 


10-1 


7.2 X 


10- 


2 


12 


42 


24, 


,0 


+32 


32 


08 


7.09 X 10- 


3 


4.78 X 


10- 


3 


9.03 X 10- 


3 


2.219 X 10" 


-2 


3.726 X 10- 


2 


6.627 X 


10-1 


1.750 X 10° 


1.36 X 


10-1 


4.0 X 


10- 


2 


12 


42 


25, 


3 


+32 


32 


50 


7.91 X 10- 


3 


5.30 X 


10- 


3 


1.106 X 10" 


-2 


2.746 X 10" 


-2 


4.000 X 10" 


2 


6.223 X 


10"! 


1.894 X 10° 


1.75 X 


10-1 


4.7 X 


10- 


2 


12 


42 


28, 


,3 


+32 


32 


44 


5.26 X 10- 


3 


3.55 X 


10- 


3 


5.52 X 10- 


3 


1.428 X 10" 


-2 


2.177 X 10- 


2 


3.789 X 


10-1 


1.189 X 10° 


1.01 X 


10-1 


3.5 X 


10- 


2 


12 


42 


30, 


,8 


+32 


32 


31 


3.32 X 10- 


3 


2.34 X 


10- 


3 


2.64 X 10- 


3 


7.98 X 10- 


3 


1.477 X 10- 


2 


2.889 X 


10-1 


9.00 X 10-1 


2.7 X 


10-^ 









o 



'^Thc 3.6 - 70 and 450 - 850 /j.m flux densities are measured in maps that have been convolved with the kernels described in Section 2.4 so as to match the higher-resolution PSFs 
to the PSF of the MIPS 160 ^m data. 

'^Calibration uncertainty: 10%; background noise uncertainty: 3 X 10~^ Jy; the "infinite" aperture corrections from Reach et al. (2005) have been applied to these data. 

'^Calibration uncertainty: 10%; background noise uncertainty: 3 X 10~^ Jy; the "infinite" aperture corrections from Reach et al. (2005) have been applied to these data. 

'^Calibration uncertainty: 10%; background noise uncertainty: 6 X 10~^ Jy; the "infinite" aperture corrections from Reach et al. (2005) have been applied to these data. 

'"'Calibration uncertainty: 10%; background noise uncertainty: 7 X 10~^ Jy; the "infinite" aperture corrections from Reach et al. (2005) have been applied to these data; note that 
these data are not stellar continuum subtracted. 

^Calibration uncertainty: 10%; background noise uncertainty: 3 X 10~^ Jy. 

^Calibration uncertainty: 20%; background noise uncertainty: 5 X 10""^ Jy. 

'^Calibration uncertainty: 20%; background noise uncertainty: 3 X 10~^ Jy. 

'Calibration uncertainty: 25%; background noise uncertainty: 1.9 X 10~^ Jy. 

^Calibration uncertainty: 10%; background noise uncertainty: 2 X 10~^ Jy; note that CO, free-free, and synchrotron emission have not been subtracted from these data, but since 
the sum of these emission components only contributes less than 10% to the global fiux density, the absence of a correction for these data should not strongly affect the results. 



-61 - 



Table 3. Global Spectral Energy Distribution of NGC 4631 



Wavelength (//m) 


Flux Density (Jy) 


Source 


Reference 


1.3 


1.56 ±0.03 


2MASS 


1 


1.6 


1.91 ±0.04 


2MASS 


1 


2.2 


1.73 ±0.04 


2MASS 


1 


3.6 


1.30 ±0.13^ 


Spitzer/IRAC 


2 


4.5 


0.83 ± 0.08^ 


Ci -I /TT> a 

Spitzer/IRAL 


2 


5.7 


2.5 ±0.3^ 


Spitzer/IRAC 


2 


8.0 


5.8 ±0.6^ 


Spitzer/IRAC 


2 


12 


5.5 ±0.8 


IRAS 


3 


24 


8.0 ±0.8 


Spitzer/MIPS 


2 


60 


83± 12 


IRAS 


3 


70 


100 ± 20 


Spitzer/MIPS 


2 


100 


210 ±30 


IRAS 


3 


160 


270 ± 50 


Spitzer/MIPS 


2 


450 


24 ±6 


JCMT/SCUBA 


(this paper) 


850 


4.5 ±0.5*^ 


JCMT/SCUBA 


(this paper) 


1230 


2.0 ±0.2*^ 


IRAM 


4 



*The "infinite" aperture corrections from Reach et al. (2005) have been 
applied to these data. 

^CO, free-free, and synchrotron emission have been subtracted from 
these global flux density measurements using data from Dumke et al. 
(2004). 



References. — 1 Jarrett et al. (2003), 2 Dale et al. (2005), 3 Rice et 
al. (1988), 4 Dumke et al. (2004) 
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Table 4. 



Results of Fitting Blackbodies Modified with A ^ Emissivity Laws to the 70 
450 fj,m SEDs of 40" Regions in NGC 4631^ 



Position (J2000) 
R.A. Dec. 



Temperature 
(K)» 



Excess 850 fim 
Flux Density (Jy) 



Ratio Excess to Total 
850 jum Flux Density 



Total Infrared 
Luminosity (Lq) 



12 


41 


40, 


,9 


+32 


31 


50 


22 


,2 


0, 


,031 ± 


0, 


,005 


0, 


,65 ± 


0, 


13 


1,33 ± 0,18, X 10* 


12 


41 


43. 


2 


+32 


31 


56 


22, 


,3 


0, 


,037 ± 


0, 


,006 


0, 


,68 ± 


0, 


,13 


1.45 ± 0.19 X 10* 


12 


41 


45. 


,8 


+32 


32 


00 


22, 


,8 


0, 


,049 ± 


0, 


,007 


0, 


,69 ± 


0, 


,13 


2.2 ± 0.3. X 10* 


12 


41 


49. 


,4 


+32 


32 


04 


23, 


,0 


0, 


,069 ± 


0, 


,009 


0, 


,65 ± 


0, 


,12 


3.4 ± 5. X 10* 


12 


41 


51. 


,8 


+32 


32 


10 


23, 


,0 


0, 


,062 ± 


0, 


,010 


0, 


,65 ± 


0, 


,12 


3.5 ± 0.5 X 10* 


12 


41 


54. 


2 


+32 


32 


14 


22, 


,6 


0, 


,058 ± 


0, 


Oil 


0, 


,56 ± 


0, 


12 


4.3 ± 0.6 X 10* 


12 


41 


,58, 


,4 


+32 


32 


12 


21 


,9 


0, 


,041 ± 


0, 


010 


0, 


,43 ± 


0, 


11 


4.2 ± 0.6 X 10* 


12 


42 


01, 


,3 


+32 


32 


17 


22 


,2 





,042 ± 


0, 


012 


0, 


,35 ± 


0, 


11 


7.0 ± 0.9 X 10* 


12 


42 


04, 


,3 


+32 


32 


25 


22 


,8 


0, 


.056 ± 


0, 


021 


0, 


,26 ± 


0, 


10 


1.6 ± 0.2 X lO" 


12 


42 


07, 


,9 


+32 


32 


35 


22 


,9 


0, 


,044 ± 


0, 


022 


0, 


,20 ± 


0, 


10 


1.9 ± 0.2 X lO" 


12 


42 


10, 


,9 


+32 


32 


35 


22, 


,6 


0, 


,034 ± 


0, 


,020 


0, 


,17 ± 


0, 


10 


1.7 ± 0.2 X 10^ 


12 


42 


14, 


,0 


+32 


32 


38 


22, 


,1 


0, 


,024 ± 


0, 


Oil 


0, 


,22 ± 


0, 


,10 


7.0 ± 0.9 X 10* 


12 


42 


16, 


,7 


+32 


32 


39 


21, 


,3 


0, 


,020 ± 


0, 


,008 


0, 


,27 ± 


0, 


,11 


3.7 ±0.5 X 10* 


12 


42 


19, 


,7 


+32 


32 


37 


21, 


,7 


0, 


,029 ± 


0, 


,007 


0, 


,41 ± 


0, 


,11 


3.1 ± 0.4 X 10* 


12 


42 


20, 


,2 


+32 


31 


59 


23, 


,0 


0, 


,015 ± 


0, 


,004 


0, 


,47 ± 


0, 


13 


1.7 ± 0.2 X 10* 


12 


42 


22, 


,0 


+32 


32 


44 


22 


,3 


0, 


,032 ± 


0, 


008 





,44 ± 


0, 


11 


3.5 ± 0.5 X 10* 


12 


42 


24, 


,0 


+32 


32 


08 


23, 


,2 


0, 


.023 ± 


0, 


005 





,57 ± 


0, 


13 


1.8 ± 0.2 X 10* 


12 


42 


25, 


,3 


+32 


32 


50 


22, 


,3 


0, 


,026 ± 


0, 


,005 


0, 


,55 ± 


0, 


13 


1.9 ± 0.3 X 10* 


12 


42 


28, 


,3 


+32 


32 


44 


22, 


,3 


0, 


.022 ± 


0, 


,004 


0, 


,63 ± 


0, 


,14 


1.17 ± 0.16 X 10* 



^The flux densities for these regions are listed in Table 2. 
'^Uncertainties in the temperatures are 2 K. 
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Table 5. Results of Fitting Semi-Empirical Models to the 5.7 - 450 /xm SEDs of 40' 

Regions in NGC 4631^ 





Position (J2000) 
R.A. Dec. 




a'' 


Excess 850 fini 
Flux Density (Jy) 


Ratio Excess to Total 
850 fim Flux Density 


Total Infrared 
Luminosity (Lq 


12 


41 40.9 


+32 


31 


50 


2.31 


0.031 ± 0.005 


0.65 + 0.13 


1.33 + 0.18 X 10' 


12 


41 43.2 


+32 


31 


56 


2.44 


0.035 ± 0.006 


0.64 + 0.13 


1.45 + 0.19 X 10' 


12 


41 45.8 


+32 


32 


00 


2.38 


0.042 ± 0.007 


0.60 + 0.12 


2.2 + 0.3 X 10* 


12 


41 49.4 


+32 


32 


04 


2.31 


0.049 ± 0.009 


0.53 + 0.12 


3.4 + 0.5 X 10* 


12 


41 51.8 


+32 


32 


10 


2.38 


0.051 + 0.010 


0.54 + 0.12 


4.3 + 0.6 X 10* 


12 


41 54.2 


+32 


32 


14 


2.31 


0.053 + 0.011 


0.50 + 0.11 


4.2 + 0.6 X 10* 


12 


41 58.4 


+32 


32 


12 


2.44 


0.039 ± 0.010 


0.41 ± 0.11 


4.2 + 0.6 X 10* 


12 


42 01.3 


+32 


32 


17 


2.31 


0.037 + 0.012 


0.30 + 0.11 


7.0 + 0.9 X 10* 


12 


42 04.3 


+32 


32 


25 


2.19 


0.039 ± 0.021 


0.18 + 0.10 


1.6 + 0.2 X 10" 


12 


42 07.9 


+32 


32 


35 


2.13 


0.030 ± 0.022 


0.14 + 0.10 


1.9 + 0.2 X lO'' 


12 


42 10.9 


+32 


32 


35 


2.25 


0.008 ± 0.020 


0.04 + 0.10 


1.7 + 0.2 X 10" 


12 


42 14.0 


+32 


32 


38 


2.31 


0.024 + 0.011 


0.22 ± 0.10 


7.0 + 0.9 X 10* 


12 


42 16.7 


+32 


32 


39 


2.50 


0.024 + 0.008 


0.31 ± 0.11 


3.7 + 0.5 X 10* 


12 


42 19.7 


+32 


32 


37 


2.44 


0.029 + 0.007 


0.41 ± 0.11 


3.1 + 0.4 X 10* 


12 


42 20.2 


+32 


31 


59 


2.50 


0.008 ± 0.004 


0.26 + 0.12 


1.7 + 0.2 X 10* 


12 


42 22.0 


+32 


32 


44 


2.25 


0.032 ± 0.008 


0.44 ±0.11 


3.5 + 0.5 X 10* 


12 


42 24.0 


+32 


32 


08 


2.50 


0.014 ± 0.005 


0.37 + 0.12 


1.8 + 0.2 X 10* 


12 


42 25.3 


+32 


32 


50 


2.56 


0.019 + 0.005 


0.41 + 0.12 


1.9 + 0.3 X 10* 


12 


42 28.3 


+32 


32 


44 


2.81 


0.016 ± 0.004 


0.47 + 0.13 


1.17 + 0.16 X 10- 



l8 



''The flux densities for these regions are listed in Table 2. 

''This is the index of the power law that describes the range of radiation flelds heating the dust. 
Uncertainties in the values are ~ 0.10. 
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Table 6. Parameters of Physical Models fit to the 5.7 - 850 jim SEDs of 40" Regions in 

NGC 4631^ 



Position (J2000) 




Total 


PAH Mass 


7^^ 


Minimum Illuminating 


R.A. 


Dec. 




Mass (Mq) 


Fraction'' 




Radiation Field'' 


12 41 40.9 


+32 


31 


50 


1.1 


X 


lO" 


0.019 


0.0151 


0.7 


12 41 43.2 


+32 


31 


56 


9.8 


X 


10^ 


0.019 


0.0103 


1.0 


12 41 45.8 


+32 


32 


00 


1.2 


X 


10" 


0.026 


0.0170 


1.0 


12 41 49.4 


+32 


32 


04 


1.3 


X 


10" 


0.026 


0.0164 


1.5 


12 41 51.8 


+32 


32 


10 


1.5 


X 


10" 


0.026 


0.0140 


1.5 


12 41 54.2 


+32 


32 


14 


1.9 


X 


10" 


0.026 


0.0153 


1.5 


12 41 58.4 


+32 


32 


12 


1.9 


X 


10" 


0.034 


0.0107 


1.5 


12 42 01.3 


+32 


32 


17 


2.7 


X 


10" 


0.042 


0.0212 


1.5 


12 42 04.3 


+32 


32 


25 


4.5 


X 


10" 


0.042 


0.0299 


2.0 


12 42 07.9 


+32 


32 


35 


3.8 


X 


10" 


0.042 


0.0278 


3.0 


12 42 10.9 


+32 


32 


35 


4.3 


X 


10" 


0.042 


0.0323 


2.0 


12 42 14.0 


+32 


32 


38 


2.2 


X 


10" 


0.042 


0.0158 


2.0 


12 42 16.7 


+32 


32 


39 


1.7 


X 


10" 


0.034 


0.0115 


1.5 


12 42 19.7 


+32 


32 


37 


1.7 


X 


10" 


0.034 


0.0172 


1.0 


12 42 20.2 


+32 


31 


59 


5.5 


X 


10^ 


0.026 


0.0110 


2.0 


12 42 22.0 


+32 


32 


44 


1.4 


X 


10" 


0.026 


0.0193 


1.5 


12 42 24.0 


+32 


32 


08 


6.2 


X 


10^5 


0.019 


0.00921 


2.0 


12 42 25.3 


+32 


32 


50 


8.1 


X 


10^ 


0.026 


0.0110 


1.5 


12 42 28.3 


+32 


32 


44 


6.6 


X 


10^ 


0.019 


0.0115 


1.0 



''The flux densities for these regions are listed in Table 2. 

''This is the fraction of the total dust mass that is in PAH molecules that are smaller than 
10^ carbon atoms. The value for the Milky way is 0.0488 ~ 0.0078. 

"^This is the fraction of dust heated by a range of radiation fields described by a power law 
distribution. The rest of the dust is assumed to be heated by a low level radiation field. 

''This is expressed in units of the local interstellar radiation field. 



Table 7. Results of Fitting Physical Models to the 5.7 - 850 /im SEDs of 40" Regions in NGC 4631^ 





Position (J2000) 




Excess 450 fim 


Ratio Excess to Total 


Excess 850 fim 


Ratio Excess to Total 


Total Infrared 




R.A. 


Dec. 




Flux Density (Jy)'' 


450 fj,m Flux Density'' 


Flux Density (Jy)*" 


850 Flux Density'^ 


Luminosity (Lq) 


12 


41 40.9 


1 QO 


Q1 

oi 


50 


— U.Uy ± U.Uo 


—0.6 + 0.3 


U.UUo ± U.UUo 


0.11 i 0.11 


1 QQ _l_ n 1 o \y 1 n8 
l.oo ±0.18 X 10 


12 


41 43.2 




Oi 


56 


— U.iU ± U.Uo 


—0.6 + 0.3 


U.Uil ± U.UUb 


n o 1 -L- nil 


1 A ^ _l_ n 1 A \y 1 n8 
1.40 it O.iy X 10 


12 


41 45.8 




oz 


00 


-U.15±U.U5 


—0.8 + 0.3 


U.UzU ± U.UU7 


A OQ _1_ AIA 
0.2<5 ± 0.10 


O O _1_ A Q 1 a8 
2.2 ± O.O X 10 


12 


41 49.4 


1 QO 


oz 


04 


— U.lo ± U.Uo 


—0.7 ± 0.3 


U.U27 ± U.UUy 


U.oU ± U.IU 


O yl 1 A C V* 1 a8 

3.4 ± 0.5 X 10 


12 


41 51.8 


+32 


32 


10 


— U.zU ± U.U7 


—0.7 + 0.3 


A r\n A _i_ n aia 
U.Uz4 ± U.UiU 


A _1_ AIA 
U.ZO ± U.iU 


A Q _1_ A 1 a8 

4.3 ± O.O X 10 


1 o 

LZ 


41 54.2 


+32 


32 


1 A 
14 


n 1 7 j_ fi 11 
— U.l ( it U.il 


— U.4 it U.o 


A A1 C _l_ A A1 1 
U.UlO it U.Uil 


A 1 i 1 AIA 
U.i4 it U.iU 


4.Z it U.o X lU 


12 


41 58.4 


+32 


32 


12 


— 0.U7 ± U.13 


—0.1 ± 0.3 


U.0U6 ± O.UlO 


A n/^ 1 n "1 n 

O.Od ± 0.10 


ylA 1 Ar'>, 1a8 

4.2 + O.D X 10 


12 


42 01.3 


+32 


32 


17 


-0.09 + 0.18 


-0.1 + 0.3 


-0.003 + 0.012 


-0.02 + 0.10 


7.0 + 0.9 X 10* 


12 


42 04.3 


+32 


32 


25 


-0.1 + 0.4 


-0.1 + 0.3 


-0.01 + 0.02 


-0.04 + 0.10 


1.6 + 0.2 X 10^ 


12 


42 07.9 


+32 


32 


35 


0.0 + 0.4 


0.0 + 0.3 


0.01 + 0.02 


0.03 + 0.10 


1.9 + 0.2 X lO** 


12 


42 10.9 


+32 


32 


35 


0.0 + 0.3 


0.0 + 0.3 


-0.01 + 0.02 


-0.07 + 0.10 


1.7 + 0.2 X 10^ 


12 


42 14.0 


+32 


32 


38 


0.1 + 0.2 


0.1 + 0.3 


-0.001 + 0.011 


-0.01 + 0.10 


7.0 + 0.9 X 10* 


12 


42 16.7 


+32 


32 


39 


0.04 + 0.14 


0.1 + 0.3 


-0.002 + 0.008 


-0.03 + 0.10 


3.7 + 0.5 X 10** 


12 


42 19.7 


+32 


32 


37 


-0.07 + 0.10 


-0.2 + 0.3 


-0.003 + 0.007 


-0.04 + 0.10 


3.1 + 0.4 X 10** 


12 


42 20.2 


+32 


31 


59 


-0.05 + 0.04 


-0.4 + 0.3 


0.004 ± 0.004 


0.12 + 0.12 


1.7 + 0.2 X 10* 


12 


42 22.0 


+32 


32 


44 


-0.08 + 0.10 


-0.2 + 0.3 


0.003 ± 0.008 


0.04 + 0.10 


3.5 + 0.5 X 10* 


12 


42 24.0 


+32 


32 


08 


-0.08 ± 0.04 


-0.6 + 0.3 


0.008 ± 0.005 


0.20 + 0.11 


1.8 + 0.2 X 10* 


12 


42 25.3 


+32 


32 


50 


-0.08 + 0.05 


-0.4 + 0.3 


0.009 ± 0.005 


0.19 + 0.11 


1.9 + 0.3 X 10* 


12 


42 28.3 


+32 


32 


44 


-0.08 ± 0.03 


-0.8 + 0.4 


0.006 ± 0.004 


0.17 + 0.12 


1.17 + 0.16 X 10* 



''The flux densities for these regions are listed in Table 2. 

''Negative numbers correspond to regions where the model fits predict a higher flux density than the observations. 



